AD-754  517 


INSTRUMENT  LANDING  SYSTEM  SCATTERING 
Chin,  L.  Jordan,  et  al 


Transportation  Systems  Center 


Prepared  for: 

Federal  Aviation  Administration 
December  1972 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Road,  Springfield  Va.  22151 


Best 

Available 

Copy 


AD 7545  1 7 


REPORT  NO.  FAA-RD-72-1 37 


INSTRUMENT  LANDING  SYSTEM  SCATTERING 


G .  Chin 
L.  Jordan 
D.  Kahn 
S.  Mor*n 

Transportation  Systems  Center 
Kendall  Square 
Cambridge,  Ma.  02142 


DECEMBER  1972 
FINAL  REPORT 


OOCUMCKT  IS  AVAILABLE  to  tKC  BUBUC 

TMBOUOM  THE  VATIOBAL  TtCMXtCM.  t 

ixromt«TioN  service.  sbbirgfielo.  .  ,  _  _  ./ 

VIRGINIA  Z2IS1. 


Prepared  for 

DEPARTMENT  OF  TRANSPORTATION 
FEDERAL  AVIATION  ADMINISTRATION 
Systems  Research  and  Development  Service 
Washington  D.C.  20591 


\!  SERVICE 


NOTICE 


This  document  is  disseminated  under  the  sponsorship 
of  the  Department  of  Transportation  in  the  interest 
of  information  exchange.  The  United  States  Govern¬ 
ment  assumes  no  liability  for  its  contents  or  use 
thereof. 


I 


I-  M*. 

FAA-RD-72-137 

•.  T.'I.  S.k'.*i. 


1  2.  Ct>r 
i 

J _ 


K*. 


INSTRUMENT  LANDING  SYSTEM  SCATTERING 


2 

G.  Chin.  L.  Jordan.  D.  Kahn  ,  S.  Morin 

9_  0'|«n>t(*>**  N«m  aM 

Department  of  Transportation 
Transportation  Svstecs  Center 
Kendall  Square,  Cambridge,  HA.  02142 

12  S»eni»Hi|  Ayec y  N«m«  *<«»* 

Department  of  Transportation 
Federal  Aviation  Administration 
Systems  Research  and  Development  Service 
Washington,  d.C.  20591 

IS  N««tt 


3,  2+t,f  mma~ i  Ci*r«|  Ha 


5„  aw»o** 

December  1972 


i  PWs»»»|  0*|^  |V  »■  ft.,,*,  N* 

D0T-TSC-FAA- 72-28 


13  »*-.  a—,  k. 

R  2103 


11  •>  N« 

FA2Q7 _ 


n  Tr„  . i  iw 
Final  Report 


1*  l|ft|  C»*r 


16  Ah*rte< 

\he  construction  of  a  mathematical  node!  of  the  Instrument 
Landing  System  (ILS)  multipath  problem  has-been  undertaken.  This 
report  presents  the  theoretical  basis  for  such  a  model,  and  newly- 
achieved  developments  in  ILS  model  construction. 


7.  k«t ILS ,  scattering 
theory,  current  deviation 
indication,  derogation,  receiver 
model,  Doppler  Shift,  DDM 


It.  Dflhb«M  S*.H. 

Availability  !S  Unlimited.  Document  may  be 
Released  to  the  Nations.'  Technical  ir.forma- 
tion  Service.  Springfield.  Virginia  22151.  for 
sale  to  the  Public. 


19.  $«Cvnty  Jttlif.  («f  Ait  report) 

20.  Sttwiff  Cltttif.  (*f  Mfl 

21-  N*.  •#  Pmpet 

22.  Pnw 

Unclassified 

Unclassified 

14® 

PREFACE 

This  final  report  of  fiscal  year  1972  on  ILS  scattering  con- 
:ains  the  new  work  performed  since  the  publication  of  the  February 
L972  report,  titled  "The  ILS  Scattering  Problem  and  Signal  Detec¬ 


tion  Model' 


In  order  that  this  reoort  be  as  self-contained  as 


possible,  most  of  the  original  work  contained  in  the  February  1972 
report  is  reproduced  here. 

Work  performed  since  the  February  report  is  found  in  Appendix 
A  which  contains  the  development  of  the  Fresnel  approximation  for 
use  in  considering  scattering  from  tall  structures;  in  Section  4 
where  the  closed-form  solution  for  scattering  from  slanted  rec¬ 
tangular  walls  is  derived  and  where  an  extension  cf  the  electro¬ 
magnetic  scattering  problem  to  non-perfectly  conducting  slabs  is 
developed;  and  in  Appendix  E  where  the  computer  program  is  dis¬ 
cussed  . 

Since  publishing  the  February  1972  report,  the  model  has 
undergone  a  validation  test  and  agrees  well  with  the  flight  data. 
The  model  has  also  been  applied  to  the  planned  Dallas  Fcrt  rtorth 
Regional  Airport  to  predict  the  expected  course  deviation  indica¬ 
tion  (CDI)  using  two  different  localizers.  These  results  may  be 
found  in  TSC  Report  Ho.  DOT-TSC-FAA-72-15,  titled  “ILS  Localizer 
Performance  Study;  Dallas  Fort  Worth  Regional  Airport  and  Model 
Validation  -  Syracuse  Hancock  Airport." 

we  gratefully  acknowledge  the  help  received  from  Mr.  S-  Lam 
and  Mr.  R.  Silva  who  contributed  their  scientific  talents  to  the 
success  of  this  research  effort,  we  would  also  like  to  thank 
Mr.  D.  Hewsom  for  his  excellent  programming  work  on  this  project. 
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Different  angles  of  alpha,  0.0*.  3.1*  and  14.0*  are 
ssed  to  position  the  specular  reflection  at  dif¬ 
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1.0  iKTROByETaGN 


Is  first  part  of  this  report ,  tbs  basic  XLS  scattering 
problem  is  investigates  starting  from  first  principles-  Harwell’s 
equations  are  first  formally  integrates  using  the  vector  Green’s 
tbsor ei.  The  resalting  integral  equations  express  tia  electric 
are  magnetic  fields  at  an  observation  point  ?  inside  a  v plane  V  in 
terns  of  volume  integrals  over  the  charge  and  carrent  distributions 
inside  ¥  j~-  sarface  integrals  of  the  electric  and  magnetic  fields 
over  the  serf  aces  bounding  ¥.  The  sarface  integrals  represent  con- 
tribctic ns  to  the  electrtstagnetic  fields  at  ?  from  radiation 
sources  locate  outside  of  the  volume  ¥  and  are  identically  zero 


sacs  ex¬ 


scar  ces. 


Tbese  integral  eqaations  are  then  applied  to  the  general 
problem  of  electromagnetic  scattering .  The  voicne  of  integration 
V  is  nade  multiply-connected  with  ta  interior  and  an  exterior 
boundary.  The  source  of  the  incident  radiation  Ce.c. ,  an  IIS 
antenna#  is  located  inside  V,  while  the  scatterer  is  enclosed  by 
the  interior  boundary  of  v.  Conseguen  tly ,  the  charge  and  current 
distributions  which  constitute  the  source  of  the  incident  radiation 
are  the  only  radiation  sources  interior  to  V.  Under  the  assumption 
that  the  perturbations  in  the  current  and  charge  distributions  of 
the  primary  source  due  to  the  presence  of  the  scatterer  can  be 


neciected 


rox- 


tsa  mi 


rations 


electric 


and  magnetic  fields  at  an  interior  point  of  V  are  applied.  The 
electromagnetic  fields  are  represented  as  suns  of  the  incident 
fields  produced  by  the  primary  source  and  the  scattered  fields 
produced  by  the  induced  currents  and  charges  in  the  scatterer-  It 
is  shown  th»t  the  scattered  electric  aid  magnetic  fields  at  the 
observation  point  P  can  be  represented  as  surface  integrals  of  the 
scattered  fields  over  the  surface  of  the  scatterer. 


To  obtain  approximate  solutions  to  these  surface  integral 
equations  for  the  scattered  fields  at  ?,  an  iterative  approach  is 
adopted .  Specifically,  frees  a  knowledge  of  the  boundary  conditions 
which  must  be  satisfied  at  the  sarface  of  the  scatterer,  approximate 
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functional  relationships  among  the  scattered  fields  and  the  known 
incident  fields  are  developed  and  then  substituted  into  the  surface 
integral  equations.  The  integrals,  which  are  now  expressed  in 
terns  of  the  known  incident  fields,  are  then  evaluated  to  produce 
approximate  expressions  for  the  scattered  fields  at  the  observation 
point.  The  functional  relationships  among  the  scattered  and  inci¬ 
dent  fields  at  the  surface  of  the  scatterer  are  extremely  compli¬ 
cated  in  the  case  of  certain  structures,  for  example ,  for  hollow 
dielectric  buildings  with  various  internal  structure,  but  very 
simple  in  the  case  of  perfect  conductors  or  buildings  with  metal 
walls  (or,  to  a  good  approximation,  metal  rod  reinforced  concrete 
vails) .  Because  of  the  difficulties  of  pursuing  the  analysis 
further  with  ncn-metal  structures  and  because  typical  airport 
scattering  objects  are  well  approximated  by  perfect  conductors,  we 
assume,  except  in  Section  IV,  that  the  scattering  objects  are 
perfect  conductors  (Approx.  No-  2). 

Application  of  the  boundary  conditions  for  perfect  conductors 
yields  a  relationship  between  the  scattered  magnetic  field  at  the 
observation  point  and  the  surface  integral  over  the  scatterer  of 
the  tangential  component  of  the  total  (incident  plus  scattered) 
magnetic  field.  To  approximate  the  total  magnetic  field  on  the 
surface  of  the  conducting  scatterer,  we  first  employ  the  principles 
of  ray  optics.  Specifically,  we  assume  as  a  first  approximation 
that  the  total  magnetic  field  is  zero  on  the  side  of  the  scatterer 
not  directly  illuminated  by  the  primary  source  (Approx.  No.  3). 

This  is  a  good  approximation  when  diffraction  effects  may  be  con¬ 
sidered  as  second  order  effects.  Diffraction  effects  may  safely 
is  considered  second  order  when  the  wavelength  of  the  incident 
radiation  is  small  compared  with  the  dimensions  of  the  scatterer. 
This  is  the  case  for  scattering  from  hangers;  however,  it  is  not 
the  case  for  scattering  from  aircraft  where  the  localizer  wave¬ 
length  and  fuselage  radius  are  comparable .  To  treat  this  case, 
special  care  would  have  to  be  taken  to  check  that  diffraction  re¬ 
mains  small;  for,  if  it  does  not.  Approximation  H^ber  3  could  not 
be  made  and  an  alternate  method  for  relating  incident  and  scattered 
fields  would  be  needed. 
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Having  assumed  that  the  tangential  component  of  the  magnetic 
field,  is  zero  or.  the  unilluminated  portion  of  the  scatterer,  it  is 
next  necessary  to  specify  it  on  the  illuminated  side.  This  is 
done  by  assuming  plane  wave  reflection  (Approx.  No.  4).  For  dis¬ 
tances  generally  encountered  in  the  ILS  problem,  this  approximation 
is  valid  (unless  the  scatterer  dimensions  are  comparable  to  the 
wavelength,  in  which  case.  Approximations  3  and  4  will  have  to  be 
modified) . 

Since  we  are  interested  in  the  values  of  the  scattered  fields 
in  the  far  field  of  the  scatterer  (the  approaching  aircraft  being 
between  the  outer  marker  and  the  far  end  of  the  runway) ,  the  in¬ 
tegral  equations  for  the  fields  may  be  expanded  asymptotically  for 
large  values  of  the  distance  between  scatterer  and  observer  (Approx. 
No.  5);  a  similar  far  field  approximation  is  made  for  tne  antenna- 
to-scatterer  distance.  Both  the  Fraunhofer  and  Fresnel  versions 
of  this  approximation  are  used  in  this  report. 

The  application  of  the  above  approximations  in  the  analysis 
leads  to  the  final  expressions  shown  in  Section  2.4,  APPROXIMATE 
SOLUTIONS  TO  THE  INTEGRAL  EQUATIONS  OF  SCATTERING  THEORY,  Equations 
(2.47)  and  (2.48)  for  che  scattered  electromagnetic  field.  These 
differ  from  those  used  by  IBM  but  are  basically  the  same  as  the 
Ohio  University  expressions.  The  differences  between  the  theory 
presented  here  and  IBM's  are  discussed  in  Sections  2.5  and  2.6, 
LOCALIZER  SIGNAL  SCATTERING  and  COMPARISON  WITH  PREVIOUS  WORK, 
where  the  theory  is  carefully  applied  to  localizer  signal  scatter¬ 
ing  by  a  rectangular  wall.  By  means  of  this  application  of  the 
theory  to  scattering  from  a  rectangular  wall,  it  is  shown  our  re¬ 
sults  reduce  to  Ohio's,  if  certain  additional  approximations  re¬ 
lating  to  reflections  from  the  ground  plane  are  made  in  our  equa¬ 
tions.  The  differences  between  our  formulation  of  the  scattering 
problem  and  the  IBM  formulation  are  of  a  fundamental  nature.  The 
practical  consequences  of  these  differences  are  shown  in  an  ac¬ 
companying  graph  in  which  the  example  used  by  IBM  for  scattering 
from  a  vertical  rectangular  wall  is  used  and  the  predicted  DDM's 
compared.  The  differences  may  often  be  significant.  For  very  tall 
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structures  even  larger  differences  would  be  found,  as  is  demon¬ 
strated  in  our  Fresnel  approximation  development,  found  in 
Appendix  A. 

In  addition  to  the  analysis  of  rectangular  wall  scattering, 
we  present  (Appendix  B)  the  new  scattering  formula  for  vertical 
triangles  and  show  how  to  use  these  when  the  triangles  are  ele¬ 
vated,  as  in  the  case  of  triangular  roof  structures  and  tail  sec¬ 
tions  of  aircraft.  We  also  present  (Section  4.4.)  a  closed-form 
solution  to  the  slanted  rectangular  wall  (not  previously  obtained 
in  the  IBM  formulation) ,  which  should  be  useful  for  calculating 
reflection  from  hangars  with  slanted  roofs.  New  closed-form 
solutions  for  double  reflection  between  two  vertical  walls  are 
given  in  Appendix  C. 

In  Section  3.0,  MODELING  OF  ILS  SIGNAL  DETECTION,  we  go  from 
the  scattering  problem  to  the  signal  detection  problem,  in  which 
we  try  to  understand  how  the  DDM  must  be  defined  in  the  presence 
of  multipath  and  develop  a  reasonable  model  of  the  ILS  signal 
detection  system.  The  IBM  and  Ohio  University  expressions  for  DDM 
are  inadequate  for  strong  multipath  environments,  since  their  ex¬ 
pressions  are  strictly  valid  only  for  single  carrier  signals  when 
the  relative  phases  of  the  received  carrier  and  sidebands  are  the 
same.  We,  therefore,  present  a  unified  model  of  ILS  signal  recep¬ 
tion  which  includes  the  dephasi.ig  of  carrier  and  sideband  signals, 
Doppler  effects,  different  receiving  antenna  gain  patterns  and 
capture  effect  systems. 

To  do  this,  a  general  expression  for  the  receiver  input  cur¬ 
rent  is  written  which  includes  the  polarization  and  gain  vector  of 
the  antenna  (the  gain  vector  for  a  small  circular  loop  receiving 
antenna  is  derived  in  Appendix  D) .  This  is  used  to  help  represent 
the  amplified  signal,  which  appears  at  the  output  of  the  IF  stage 
of  the  receiver,  in  terms  of  the  different  transmitted  modulation 
waveforms  and  the  gains  and  phase  delays  associated  with  the  dif¬ 
ferent  radiation  paths.  This  IF  signal  is  then  passed  to  a  second 
detector  which  generates  an  audio  frequency  signal  which  is  passed 
through  a  set  of  filters  to  obtain  the  relative  90  Hz  and  150  Hz 
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amplitudes,  from  which  the  course  deviation  indication  (CDI)  may 
be  determined. 

Since  the  second  detector  (amplitude  modulation  detector)  is 
a  nonlinear  device,  its  output  reflects  interactions  between  the 
intended  ILS  signal  and  the  spurious  signals  received  due  to  multi- 
pathing  or  due  to  transmission  of  a  secondary  carrier.  In  order 
to  estimate  the  relative  passage  of  this  output  through  the  selec¬ 
tive  SG  Hz  and  150  Hz  tone  filters,  Fourier  analysis  is  used  to 
express  the  detected  audio  signal  in  terms  of  discrete  frequency 
components  for  which  well  defined  transmissivities  by  each  filter 
are  assumed.  The  interaction  between  the  course  and  clearance  sig¬ 
nals  in  the  general  case  of  a  dual  carrier  frequency  system  gives 
rise  to  the  much-utilized  "capture  effect".  Advantage  is  taken  of 
the  large  separation  between  the  ILS  signal  modulation  frequencies 
and  the  intercarrier  beat  (8kHz)  to  find  an  approximate  linear  ex¬ 
pression  for  the  total  detected  audio  signal  in  terms  of  independ¬ 
ent  course  and  clearance  audio  signals. 

The  isolated  course  and  clearance  s. ;nals  have  the  character 
of  audio  outputs  from  an  AM  detector  generated  by  standard  single 
carrier  ILS  signals,  distorted  by  multipathing .  For  this  single 
carrier  case,  a  simple  relation  is  found  for  the  principal  com¬ 
ponents  of  the  detected  audio  signal  lying  within  the  passbands  of 
the  modulation  frequency  filters.  The  analysis  in  this  case  is 
valid  in  the  approximation  that  the  aggregate  of  interfering  sig¬ 
nals  is  somewhat  weaker  than  the  direct  carrier  signal  -  a  reason¬ 
able  condition  for  any  marginally  flyable  course. 

Doppler  modifications  of  detected  signals  arise  as  a  conse¬ 
quence  of  the  variation  of  multipath  phase  delays  which  is  implicit 
with  the  motion  of  an  airborne  receiver.  In  the  approximation 
stated  immediately  above,  each  received  component  of  multipath 
interference  may  be  characterized  by  the  relative  Doppler  shift  of 
its  carrier  to  the  direct  path  carrier.  Sum  and  difference  com¬ 
binations  of  modulation  and  Doppler  frequencies  can  under  certain 
circumstances  give  net  frequencies  near  the  filter  center  fre¬ 
quencies  resulting  in  possible  false  signals  being  passed  by  the 
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90  Hz  or  150  Hz  filters.  With  the  aircraft  approach  speeds  used 
today,  the  Doppler  may,  for  example,  increase  the  90  Hz  signal  to 
the  point  where  it  passes  through  the  150  Hz  filter,  or  it  may  in¬ 
crease  the  150  Hz  to  the  extent  that  it  is  excluded  from  the  150  Hz 
filter.  These  possibilities  are  investigated  by  studying  the  fre¬ 
quency  response  of  a  narrowband  modulation  filter  and  calculating 
all  significant  contributions  in  each  filter  output.  Thus,  values 
are  obtained  for  the  two  detected  modulation  amplitudes  from  which 
the  apparent  DDM  and  CD1  are  conventionally  determined.  As  a 
practical  matter,  it  appears  that  the  Doppler  contributions  to  sig¬ 
nal  derogation  cannot  be  neglected  for  aircraft  speeds  in  the  vicin¬ 
ity  of  200  feet/sec.  and  filter  bandwidths  around  15  cps. 

Thus  the  new  model  is  capable  of  calculating  the  detected  audio 
frequency  signal,  including  Doppler  and  capture  effects,  in  the 
presence  of  any  moderately  strong  multipath  interference.  It  may 
therefore  be  applied  to  determination  of  the  apparent  CDI  for  any 
conventional  localizer  or  glide  slope  system. 

In  Section  4.0  we  present  a  formulation  of  the  electromagnetic 
scattering  problem  for  imperfectly  conducting  thin  flat  slabs  as  a 
possible  extension  of  the  theory  presented  in  Section  2.0  for  use 
in  situations  where  the  electromagnetic  waves  are  incident  upon 
non-metallic  structures  at  large  angles  of  incidence. 

Finally,  in  Appe.  dix  E,  the  computer  program  itself  is 
discussed. 
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2.0  THE  ILS  SCATTERING  PROBLEM 


2.1  INTRODUCTION 

Electromagnetic  scattering  is  one  of  the  most  important  and 
most  complex  topics  in  classical  mathematical  physics.  Basically, 
the  problem  is  to  determine  the  field  perturbations  (scattered 
fields)  which  ' jsult  when  various  obstacles  (the  scatterers)  are 
placed  in  known  (incident)  electromagnetic  fields.  To  solve 
a  given  scattering  problem,  one  attempts  to  find  a  solution  of 
Maxwell's  equations  which  has  the  property  that  when  it  is  added 
to  the  known  or  incident  electromagnetic  field  the  resulting  total 
field  satisfies  the  appropriate  boundary  conditions  at  the  surfaces 
cf  the  scattering  obstacles.  Unfortunately,  exact  solutions  to 
scattering  problems  have  been  obtained  in  only  a  limited  number  of 
cases.  In  this  report,  integral  equations  for  the  scattered  fields 
will  be  developed  by  directly  integrating  the  electromagnetic 
field  equations.  Approximate  solutions  to  these  integral  equations 
for  the  case  of  perfectly  conducting  scatterers  will  then  be 
applied  to  investigate  Jocalizer  signal  scattering  by  a  flat, 
vortical  wall.  Comparisons  will  be  made  with  previous  work. 

2.2  INTEGRATION  OF  MAXWELL'S  EQUATIONS 

The  Maxwell's  equations  in  M.K.S.  units  for  a  homogeneous, 
isotropic  medium  with  permittivity  e  and  permeability  y  are  given 
below 


7  x  E  =  iu;uH 


7  x  H  =  J-iweE 


7  *  H 


7  •  E  =  P/c 


(2.1) 

(2.2) 

(2.3) 


(2.4) 


In  Equations  (2.1)  through  (2.4),  it  has  been  assumed  that  all 
fields  vary  harmonically  in  time  as  e  The  quantities  E,  H,  J, 

and  f»  are,  respectively,  the  electric  field,  the  magnetic  field, 
the  current  distribution,  and  the  charge  distribution.  Taking  the 
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curls  of  Equations  (2.1)  and  (2.2)  we  find  that  E  and  H  satisfy 
the  following  field  equations: 

2 

7  x  (7xE)  -  k  E  =  iuuJ  ,  (2.5) 

7  x  (7xK)  -  k  H  =  7xJ  ,  (2.6) 

2  2 

where  k  =  u  eu. 

Let  V  be  a  closed  volume  in  the  medium  bounded  by  a  regular 
surface  S.  Let  Q  and  P  be  vector  fields  defined  at  each  point  in 

•f  -► 

the  medium.  The  vector  Green's  theorem  for  P  and  Q  has  the 
following  form: 

J  (Q*7x7xP  -  P*7x7xQ)dv 
V 

=  f (Px7xQ  -  QxVxP) -n  ds  ,  (2.7) 

where  n  is  the  unit  normal  to  the  surface  S  and  is  directed  out¬ 
ward  from  the  volume  V. 

The  integral  identity  given  by  Equation  (2.7)  can  be  used  to 
directly  integrate  the  field  Equations  (2.5)  and  (2.6).  To  begin 
with,  we  define  the  function  ?(r',r)  as  follows: 

ikl?*-?! 

? (r*  ,r)  =  e _  .  (2.8) 

!?•-*! 

The  two  point  function  ?  is  just  the  Green's  function  of  the 
Helmholtz  equation.  That  is  to  say,  y  satisfies  the  equation 

2  2 

(7+k  )?  =  -4:t5  (r'-r)  ,  (2.9) 

where  6  (r'-r)  is  the  Dirac  Delta  Function.  The  vector  r'  can  be 
thought  of  as  the  position  vector  relative  to  some  origin  0  of  an 
observation  point  P  inside  V  while  r'  is  the  position  vector  rela¬ 
tive  to  the  origin  of  any  source  point  in  V  or  on  S.  Following 
Stratton,  we  first  define  the  vectors  P  and  Q  appearing  in 
Equation  (2.7)  as  follows: 
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*  =  £(r) 

Q  =  .a  'UZ'.r) 


(2.10) 


where  a  is  a  unit  vector  in  some  arbitrary,  fixed  direction. 
Referring  to  Equations  (2.5),  (2.9),  and  (2.10),  the  following 
vector  relations  hold: 

$x<5  =  5?xa 

7x{7xQ)  =+ak2ij>+  4ta  a  (r*-r)  +  7(a*7?)  (2.11) 

-+■  -*>  0-*-  “•* 

7x(7xP)  =  It  E  +  is.uJ 

Substituting  these  various  expressions  into  the  integral 
relation  (2.7)  and  integrating  with  respect  to  the  unprimed  source 
coordinates,  we  obtain  after  a  few  trivial  vector  manipulations 
the  following  expression  for  E (r  * ) 


E(r’)  =  ^(ioyjy-i-payjdv 


1  f 

>  -"(ti 

~4rr/|ia 


7T/  iuy (nxH) ¥+ (nxE)x7¥+ (n*E) V?  ds 


E)V?J« 


(2.12) 


The  unit  vector  a  has  been  dropped  from  (2.12)  since  its  orienta¬ 
tion  is  arbitrary  and  it  would  simply  scalar  multiply  both  sides 
of  (2.12).  Repeating  this  analysis  with  P  defined  as  H(r),  the 

-f 

following  expression  for  H(r  )  is  obtained: 


■+  -*•  1  f  +  i 

H  (r* )  =  47  J  (Jx7¥)dv 


(2.13) 


Lf  r 

4tt  * 


(nxE) ?- (nxH) x7¥ 


-  -  .  .1 

7¥-(n*H)7?J< 


Equations  (2.12)  and  (2.13)  express  the  electric  and 
magnetic  field  intensities  at  any  point  P  inside  the  volume  V  in 
terms  of  volume  integrals  over  the  charges  and  currents  inside  V 
and  surface  integrals  of  the  fields  over  the  bounding  surface  S. 
The  surface  integrals  represent  the  contributions  to  £  and  £  from 
sources  located  outside  V.  These  two  equations  will  be  the  basis 
for  our  treatment  of  elec.  aagnetic  scattering. 
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2.3  INTEGRAL  EQUATIONS  IN  THE  THEORY  OF  ELECTROMAGNETIC  SCATTERING 


To  begin  our  discussion  of  electromagnetic  scattering,  con¬ 
sider  the  geometry  depicted  in  Figure  2.1.  Figure  2.1  depicts  a 
multiply-connected  volume  V  containing  a  localized  source  of  elec¬ 
tromagnetic  radiation  which  occupies  a  volume  in  V.  This  source 
might,  for  example,  be  a  glide  slope  or  localizer  antenna.  The 
charge  and  current  distributions  within  V.  will  be  denoted  by  p. 
and  J^,  respectively.  Elsewhere  in  V,  J  and  p  are  assumed  to  be 
identically  zero.  Note  that  the  volume  V  is  bounded  by  an  interior 
boundary  surface  S^  and  an  outer  boundary  surface  SQ.  The  volume 
V1  bounded  by  S^  is  not  included  in  V.  The  volume  V  consists  of 
the  region  bounded  by  SQ,  excluding  the  interior  of  S1 .  To  apply 
Equations  (2.12) 


and  (2.13)  to  such  a  multiply-connected  region,  the  surface  inte¬ 
grals  must  be  carried  out  over  both  bounding  surfaces. 


Let  us  assume  that  the  source  occupying  the  volume  in  V  is 

the  only  source  of  radiation.  Let  E^(r)  and  H^(r)  denote,  respec¬ 
tively,  the  electric  and  magnetic  field  intensities  produced  by 
this  source.  The  subscript  "i"  is  used  to  denote  incident  field. 
The  values  of  and  at  a  point  P  inside  V  are  given  by  the 
following  equations  (Eqs.  (2.12)  and  (2.13)). 


E±  (r*  )  =  1_  f (i*uJir+  1  c^rjdv 
4"  v  £ 

-  ^  (nxHi>;+(nxEi)xy?+(n-Ei)??J 


ds  (2.14) 


r-  S- 
3=0  3 


2,  (r* )  =  (5  x$?)  dv 

x  y  i 

l_y^  J  j^iac  (nxEii  ?-  (nxH^xVy-  (n-Hi) ds. 


j=o  Sj 


(2.15) 


Note  that  in  Equations  (2.14)  and  (2.15),  the  volume  integrals  are 
extended  only  over  the  volume  V.  occupied  by  the  localized  source 
since  J  fc’-u  o  are  zero  elsewhere  :n  V.  Note  also  that  each  equa¬ 
tion  contains  two  surface  integrals  because  of  the  two  surfaces 
which  bound  V. 

Equations  (2.14)  and  (2.15)  represent  solutions  to  the  field 
Equations  (2.5)  arid  (2.6)  for  the  electric  and  magnetic  field  dis¬ 
tributions  inside  the  multiply  connected  region  due  to  the  local¬ 
ized  source  occupying  the  volume  in  V  depicted  in  Figure  2.1. 
Since  this  localized  source  is  assumed  to  be  the  only  source  of 
radiation  in  the  universe,  tie  surface  integrals  appearing  in 
Equations  (2.14)  and  (2.15)  are  in  fact  identically  zero.  To  see 
that  this  must  be  so,  consider  first  the  surface  integrals,  in 
Equations  (2.14)  and  (2.15)  of  and  iL  over  the  interior  boundary 


S^.  Since  encloses  no  sources  (J 


0  and 


=  0  in  Vj),  we  can 
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contract  the  surface  down  to  a  new  boundary  surface  S,  which  is 
enclosed  by  without  adding  any  new  charges  or  currents  to  the 
volume  of  integration-  Consequently,  the  only  changes  in  Equa¬ 
tions  (2.14)  and  (2.15)  which  result  free  shrinking  S.  down  to  S, 
is  that  the  surface  integrals  of  E^  and  over  S.  are  now  carried 
out  over  S^.  Since  the  electric  and  magnetic  field  intensities  at 
?  do  not  change  as  a  result  of  this  contraction,  we  must  conclude 
that  the  surface  integrals  over  S,  equal  the  corresponding  surface 
integrals  over  S^.  But,  the  values  of  the  surface  integrals  in 
(2.14)  and  (2.15)  of  E^  and  1L  ever  can  he  made  arbitrarily 
small  by  simply  allowing  the  surface  area  of  to  go  to  zero. 
Consequently,  ve  can  conclude  that  the  surface  integrals  of  E^ 
and  over  appearing  in  Equations  {2.14}  and  (2.15)  are 
identically  zero.  That  is  to  say. 


and 


jf  |i<ub  (nxiL )?  +  tnxii } x7y+  (n- E- )  7?Jds  = 


(2.163 


(nxEi)?-(nxHi)x7»-(n*Hi)f?j 


ds  =  0  (2.17) 


h  similar  argument  can  be  made  regarding  the  surface  integrals 

of  E -  and  H.  over  the  outer  boundary  S  .  Since  no  new  sources  of 
1  jl  o 

radiation  are  added  to  the  volume  c-f  integration  by  expanding  s 

to  seme  new  surface  S  ( j  =  o  and  o  2  o  outside  S  } ;  the  surface 
_  o  o 

integrals  over  SQ  must  equal  the  corresponding  integrals  over  Sc- 

But  if  the  surface  SQ  is  allowed  to  recede  to  infinity,  the  surface 

integrals  ever  S  go  to  zero  because  of  the  following  asvsototic 
-  o  '  • 

orooerties  of  electromagnetic  fields  c reduced  by  localized 

1,2 


sources. 


MMl***" 


and 


J~ fi«£  0. 

5  L  J 


C2.SU 


Substituting  Equations  (2.161,  (2.17) ,  (2.20)  and  (2.21)  into 
bioas  (2.14)  and  (2.15) ,  Me  obtain  the  following  expressions 


-  >  r 

**  1 

t? 


iL  Cr*)  =  f  f5-*$rldv 


(2.22) 


(2.23) 


Equations  (2.22)  ana  (2.23)  are  the  fan: liar  expressions  for  the 
electric  and  magnetic  field  intensities  produced  by  localized 
sources  of  radiation  in  the  absence  of  boundaries. 


Suppose  now  that  a  scattering  coject  is  placed  inside  the 
interior  boundary  S, .  The  Situation  is  depicted  in  Figure  2.2. 


Figure  2.2.  Geometry  Illustrating  the  Arrangement  of 
Source  and  Scatterer 

The  scatterer  could  be  cosssosed  of  a  conducting  or  dielectric 
malarial  or  a  combination  of  the  two.  Chaxv,.^,  and  currents  vill 
be  induced  in  the  scatterer  bv  the  incident  fields  E -  and  H- 

j.  X 

generated  by  the  urinary  source  of  radiation  occupying  the  volume 
V.  in  V.  these  secondary  sources  will  in  turn  radiate  electro¬ 
magnetic  saves.  Let  Es (r)  and  Hs (r)  denote,  respect ively ,  the 
electric  and  magnetic  fields  generated  by  the  charges  and  currents 
induced  in  the  scatterer.  the  subscript  "s*  is  used  to  denote 
scattered  field,  the  total  electric  and  magnetic  fields  can  be 
represented  as  follows: 

ECr)  =  E±  Cr)  +  E. Cr) 

_  -  -  -  (2.24 

H(r)  =  H£  Cr)  *  Hs(r) 

the  scattered  fields  E  ar.d  H  will  of  course  exert  forces 
on  the  charges  and  currents  of  the  primary  source.  As  a  result. 


the  charge  distribution  and  the  current  distribution  J.  will 

assure  new  values  which  we  will  denote  by  p£  and  ,  respectively. 

However,  we  will  assure  that  the  perturbations  of  J.  and  p.  due  to 

1  x 

the  presence  of  the  scatterer  are  second  order  effects  and  that 
they  can  be  ignored.  That  is  to  say,  we  will  assure  that 


pi 


(2.25) 


These  approximations  should  be  very  good  for  scatterers  in  the  far 
field  of  the  urinary  radiator.^ 


To  calculate  the  total  electric  and  magnetic  fields  at  a 
point  P  inside  V,  we  s imply  apply  Equations  (2.12)  and  (2.13).  For 
example,  the  electric  field  is  given  by 


£(r*>  =  J~ (iasyjjv+i  -^7?)  dv 

V. 

x 


iwc(nxH}?-MnxE)x%+{n-E}f?  ds 


(2.26) 


where  E  =  E-  +  E  and  3  =  H.  +  E_  (Eq.  (2.24)).  Notice  that  the 
localized  noui ce  occupying  che  volume  V.  is  still  the  only  source 
of  radiation  inside  V  since  the  scatterer  lies  in  the  volure 
bounded  by  S1  and  this  region  is  excluded  f ror  V.  Notice  also 
that  we  have  made  use  of  the  approximation  (2.25) .  According  to 
Equation  (2.22),  the  volume  integral  appearing  in  Equation  (2.26) 
equals  E-(r’).  Consequently,  since  E  =  E^  +  Es,  the  scattered 
electric  field  at  P  is  given  by 


*  ifr 

3=o  3 


i  (nxH.  )f  +  (nxE. )  x 


5?+(n-E.  )T?J 


cs 


1_ 


3=o  Sj' 


i-u  (nxK  5?  +  (nxE  ) xv'r+  (n 

5  S 


V5*] 


ds  (2.27) 
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where  we  have  made  use  of  Equation  (2.24)  and  have  separated  the 

surface  integrals  involving  the  incident  and  scattered  fields. 

But  the  integrals  of  E-  and  H,  over  the  surfaces  S  and  S,  have 

already  been  shown  to  be  identically  zero  (Eqs.  (2.16)  and  (2.23)}. 

Furthermore ,  if  the  surface  S  is  allowed  to  recede  to  inf in it v. 

o 

the  surface  integral  over  Sq  of  the  scattered  fields  will  go  to 

zero  since  E  and  B  have  the  same  asvnototic  orooerties  as  E.  and 
„  s  s  ....  i 

H-  (Eqs.  (2.18)  and  (2.19)).  Consequently,  we  finally  arrive  at 
the  following  expression  for  (r ' ) : 


Es(r')  = 


-  IT 


(nxH  )y+(ftxE  )xV?+(n 
s  s 


(2.28) 


Since  the  surface  SQ  has  been  allowed  to  recede  to  infinity. 
Equation  (2.28)  is  valid  for  any  point  P  lying  outside  the  surface 
S. .  The  analogous  expression  for  H  (r*)  is 

x  S 

Hs  (r>)  =  57  y^jiwe(nxEs)?~(nxHs)xVf~{n*Es)?fJ  ds. 


(2.29) 


In  fact,  the  surface  can  be  collapsed  down  onto  the  surface  of 
the  scatterer  itself  so  that  Equations  (2.28)  and.  *2.29)  become 
valid  for  all  points  P  outside  the  body  of  the  scatterer.  This 
situation  is  depicted  in  Figure  2.3. 


SCATTERER 


Figure  2.3.  Collapsing  Surface 
S ^  onto  Scatterer 


Henceforth  then,  it  will  be  un¬ 
derstood  that  the  surface  S^  is 
the  actual  surface  of  the  scat¬ 
terer.  The  unit  normal  n  which 
appears  in  Equations  (2.28)  and 
(2.29)  points  into  the  body  of 
the  scatterer. 

The  integral  Equations 
(2.28)  and  (2.29)  express  the 
scattered  electromagnetic  fields 
at  an  observation  point  P  in 
terms  of  integrals  of  the 
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scattered  fields  over  the  surface  of  the  scatterer.  These  equa¬ 
tions  represent  first  order  solutions  to  the  scattering  problem 
in  that  they  are  based  upon  the  approximation  that  the  charge  and 
current  distributions  of  the  primary  source  are  not  changed  by  the 
fields  generated  by  the  scatterer  (Equation  (2.25)).  In  the  next 
section,  useful  approximate  solutions  to  Equations  (2.28)  and  (2.29) 
will  be  developed  for  the  case  of  perfectly  conducting  scatterer.?. 


2.4  APPROXIMATE  SOLUTIONS  TO  THE  INTEGRAL  EQUATIONS  OF  SCATTERING 
THEORY 

Approximate  solutions  to  the  integral  Equations  (2.28)  and 

(2.29)  can  oe  obtained  by  using  iterative  methods.  Specifically, 

from  a  knowledge  of  the  boundary  conditions  which  must  be  satisfied 

on  the  surface  of  the  scatterer  S. ,  approximate  functional  rela- 

tionships  among  Eg  and  Hs  .  n  S^  and  the  known  incident  fields  E^ 

and  on  S^  are  developed  and  then  substituted  into  the  integrals 

in  Equatiors  (2.28)  and  (2.29).  The  integrals  are  then  evaluated 

to  produce  approximate  expressions  for  E  and  H  at  an  observation 

s  s 

point  P.  This  procedure  is  particularly  straightforward  for  the 
case  of  perfectly  conducting  scatterer s  because  of  the  relative 
simplicity  of  the  boundary  conditions.  In  this  section  we  will  con¬ 
centrate  exclusively  upon  perfectly  conducting  scatterers. 

Consider  the  integral  Equation  (2.29)  for  the  scattered 
magnetic  field  intensity: 


H  ir') 
s 


(nxEs,'j)-  (nxHs)x^'t’-  (n* 


ds 


(2.29) 


For  convenience.  Equation  (2.29)  will  be  modified  by  making  use 
of  the  fact  that,  as  shown  earlier  the  surface  integral  appearing 
in  (2.29)  of  the  incident  fields  E^  and  EL  is  identically  zero 
(Eq.  (2.17)).  That  is  to  say. 


0 


1_ 

4tt 


j'  |i^£  (nxE^)F- (rxHi)x$,{'- (n 


(2.30) 
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Adding  Equation  (2.30)  to  (2.29)  and  recalling  that  E  =  E.  +  E  and 
-*■-*•*  is 

H  =  Hi  +  Hs,  we  obtain  the  following  equation: 


H 


£')  =  f  |iwe(nxE)y-(nxH)xfa-(n‘H)faJds,  (2.31) 


where  E  and  H  are  the  total  fields  on  S^.  At  the  surface  of  a 
perfect  conductor,  the  tangential  component  of  the  total  electric 
field  and  the  normal  component  of  the  total  magnetic  field  are 
identically  zero.  Mathematically,  these  boundary  conditions  can 
be  expressed  as  follows: 

n  x  E  =  0;  (nx2  =-nxE. ) 

s  r 


n  •  H  =  0;  (n*H  =-n*H. ) 

s  1 


(2.32) 


Consequently,  for  a  perfect  conductor.  Equation  (2.31)  takes  the 
following  simple  form: 


Hs(r')  =  - 


1_ 

4ir 


/ 


(nxH)  x  'fa  ds 


(2.33) 


To  approximate  (nxH)  on  the  surface  of  the  scatterer,  we  first 
employ  the  principles  of  geometrical  or  ray  optics.  According  to 
ray  optics,  there  is  a  deep  shadow  region  (no  illumination)  on  the 
side  of  the  scatterer  not  directly  exposed  to  the  incident  radia- 
tion  from  the  primary  source.  In  the  shadow  region,  E  and  H  are 
identically  zero  (E  =-E.  and  H  =-H.).  The  surface  S.  is  thus 
divided  into  an  illuminated  side  S+  and  a  shadow  side  S_,  the  two 
sides  being  separated  by  a  sharp  shadow  boundary  F .  This  situation 
is  depicted  in  Figure  2.4.  If  it  is  assumed  that  E  and  H  are 
identically  zero  on  S  ,  then  Equation  (2,33)  becomes 


Hs(r') 


=  ~  h  f 


(nxH)  x  fa  ds 


(2.34) 
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SHADOW 


SOURCE 


Figure  2.4  Shadow  and  Illuminated  Regions  of  Scatterer 

where  the  integral  is  performed  only  over  the  illuminated  side  of 
the  scatterer.  Finally,  to  specify  (nxS)  on  S+,  the  assumption  of 
plane  wave  reflection  from  the  local  infinite  tangent  plane  is 
made.  That  is  to  say,  at  each  point  Q  on  S+,  it  is  assumed  that 
E  (Q)  and  H  (Q)  have  the  values  of  the  reflected  fields  which 

-t- 

would  exist  if  a  plane  electromagnetic  wave  with  amplitudes  E. (Q) 
and  IK  (Q)  were  incident  upon  an  infinite  plane  perfectly  conducting 
boundary  tangent  to  S+  at  Q.  The  boundary  conditio:  satisfied  by 
the  tangential  components  of  the  incident  and  reflected  magnetic 
fields  when  a  plane  wave  is  incident  upon  an  infinite  ,  plane  , 
perfectly  conducting  boundary  is 


n  x  H  =  fi  x  H.  , 
s  i 


'2.35) 


or  equivalently. 


n  X  H  =  n  x  (H^+Hg)  =  2 (nxH^) 


(2.36) 


Assuming  that  Equation  (2.36)  is  valid  at  each  point  on  S+, 
Equation  (2.34)  becomes 


Hs( 


r ' )  =  -  /  (nxH. )  x  ds 

«4”  J  1 


(2.37) 


The  analogous  expression  for  Es(r')  can  be  obtained  by  apply¬ 
ing  Equation  (2.2)  which,  in  the  absence  of  currents,  takes  the 
form 
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?  x  H  =  -  it 


ltlE  E 


(2.38/ 


Substituting  (2.37)  into  (2.38)  the  following  expression  for  E^r') 
is  obtained: 

V?')  =  -  5^37  ?'  x  f  (nxHi)  x  ?y  ds  ,  (2.39) 


where  9*  acts  upon  the  primed  or  field  point  coordinates.  Equation 

•f 

(2.39)  can  be  rewritten  by  taking  the  operator  9*x  inside  the 
integral  and  noting  that  V?  is  the  only  function  under  the  integral 
which  depends  upon  r': 

=  '  2ruTe  f  [(nx2i)^'-  ($<?)  -  {(nxH^-S'}  Wlds 

s+  (2.40) 

-f-  4, 

Since  9V  is  a  function  of  |r*-r|,  the  following  relationships  hold 


•  (?y)  >  ($y)  =  -  v2y 

j^(nxSi)*9«J  -  ^{nxHj.^j  Vy 


(2.41) 


Substituting  the  relations  (2.41)  into  (2.40)  and  noting  that  for 
rVr,  -?2'?  =  k2?  (Eq.  (2.9)),  we  obtain  the  following  equation 


for  Eg (r ' ) : 


(r'>  =  ~  37  (?)  /  J  [^xHx)y+  —  |(HxHi)*:5}^y  ds 

'  fi  *  ^10  AO'S 


(2.42) 


where  we  have  made  use  of  the  fact  that  k=u/ye. 

We  will  be  primarily  concerned  with  the  distributions  of  Es 

and  H  in  the  far  field  or  Fraunhofer  zone  of  the  scatterer,  To 
s  -+ 

obtain  expressions  for  E  and  H  in  the  Fraunhofer  zone,  Equa- 

5  5 

tions  (2.37)  and  (2.42)  are  expanded  asymptotically  for  large 
values  of  r*=|r,j.  For  convenience,  let  the  origin  of  coordinates, 
0,  lie  on  the  surface  of  the  scatterer.  Let  D  denote  some  charac¬ 
teristic  linear  dimension  of  the  scatterer.  The  Green's  function 


WPIS 


*?  is  given  in  Equation  {2.8): 


eik|r*-r|  ■ 


(2.8) 


Suppose  that  r*>>D.  Then  in  the  denominator  of  (2.8),  |r'-r|  can 
be  approximated  very  accurately  by  r*.  That  is,  the  variations  in 
the  amplitude  of  1  as  r  varies  over  the  surface  of  the  scatterer 
can  be  ignored  if  r*>>D.  More  care  must  be  exercised  in  approxi¬ 
mating  the  phase  of  V  because  of  the  oscillatory  behavior  of  the 
exponential.  If  r*>>D,  then  |r'-rj  can  be  approximated  to  terms 
of  second  order  in  r  =  |r|  as  follows: 


2 

r '  -  r 1  *  r  +  - — 
2r ' 


(r'-r)2 

2r  ’  ' 


(2.43) 


where  r'=r’/r'  is  a  unit  vector  in  the  direction  of  r*.  In  the 

2 

Fraunhofer  zone,  D  /2r*  is  very  small  compared  to  the  wavelength 
>.=2^/k  of  the  incident  radiation  so  that  the  quadratic  terms  in 
Equation  (2.43)  can  be  ignored  in  comparison  with  the  other  two 
terms: 


=s  r»  -  r*-r  (Fraunhofer  approximation)  (2.44) 

Combining  these  amplitude  and  phase  approximations,  we  may  write 
down  the  following  asymptotic  expression  for  ?  for  field  points  in 
the  radiation  zone  of  the  scatterer: 


e-ik(r’-r) 


(2.44A) 


A  similar  analysis  of  V?  and  (A*V)^V  where  A  is  an  arbitrary  vector 
leads  to  the  following  asymptotic  expressions: 


$¥  =  -  ik 


e-ik(r’  *r)£, 


(2.45) 


(A-$)$* 


(A*r' ) r'k2 


ikr’ 


-ik (r' *r) 


(2.46) 


Substituting  the  asymptotic  forms  (2.44),  (2.45)  and  (2.46) 
into  Equations  (2.37)  and  (2.42),  we  obtain  the  following  asymptotic 
or  far  field  expressions  for  Hg  and  Eg: 
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ds,  (2.47) 


Hs(r’)  = 


(nxH^)  e 


-ik  (r '  *r) 


es (i ' )  -  ^(t)V2  •  r'  x  ?■*  /  (Kxiru-^'^'as, 

'  *  L  S  (2.48) 


where  we  have  made  use  of  the  identity  a  x  (axb)  =  (a-b)a  -  a  b  in 
Ecuation  (2.48).  Notice  that  the  asvmototic  excressions  for  H 

s 

and  £s  in  Equations  (2.47)  and  (2.48)  are  related  as  follows: 


r'  x  H  + 
s 


1/2  _ 

E  =  0 
s 


r’  x  E  -  H  =0 
s  s 


(2.49) 


These  results  are  in  agreement  with  the  asymptotic  relations 
given  in  Equations  (2.18)  and  (2.19). 

Equations  (2.47)  and  (2.48)  are  the  basic  equations  used  by 
3 

the  Ohio  University  researchers  and  the  TSC  group  for  studying 
localizer  signal  scattering  from  metallic  airport  structures  with 
plane  surfaces.  However,  caution  must  be  exercised  in  applying 
these  equations  since  the  assumptions  upon  which  they  are  based 
are  unrealistic  for  many  applications.  Generally  speaking. 

Equations  (2.47)  and  (2.48)  can  be  applied  with  some  confidence  in 
situations  in  which  the  wavelength  X=2-/k  of  the  incident  radiation 
is  small  compared  with  the  dimensions  of  the  scatterer.  When  X 
is  comparable  to  or  greater  than  the  dimensions  of  the  scatterer 
(length,  radius  of  curvature,  etc.)  the  assumption  of  a  sharp 
boundary  on  the  surface  of  the  scatterer  separating  the  illuminated 
side  (S^)  from  the  shadow  side  (S_)  breaks  down  because  of  dif¬ 
fraction.  The  assumption  of  local  plane  wave  reflection  from  the 
infinite  tangent  plane  also  breaks  down  when  X  is  comparable  to  or 
greater  than  the  radius  of  curvature  of  the  scatterer.  It  would 
for  example,  be  unwise  to  attempt  to  use  Equations  (2.47)  and  (2.48) 
to  study  localizer  signal  scattering  from  the  fuselages  of  aircraft 
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since  at  localizer  frequencies,  A =10  ft  which  is  comparable  tc  or 

greater  than  the  radii  of  curvature  of  all  aircraft  presently  in 

use.  Another  fact  to  be  borne  in  mind  is  that  Equations  (2.47) 

and  (2.48)  apply  to  the  Fraunhofer  zone  of  the  scatterer.  That 

is  to  say,  these  equations  are  based  upon  the  assumption  that 
2 

D  /2r*  is  very  small  compared  to  a  wavelength  where  D  is  some 
characteristic  linear  dimension  of  the  scatterer.  For  very  large 
scatterers,  and  for  observation  points  relatively  close  to  the 
scatterer,  this  condition  may  not  be  satisfied.  One  way  of  cir¬ 
cumventing  this  difficulty  is  to  divide  up  the  illuminated  surface 
S+  into  smaller  subsections,  apply  Equations  (2.47)  and  (2.48)  to 
each  subsection,  and  then  add  up  the  contributions  from  the 
various  subsections  to  get  the  total  scattered  fields.  A  more 
direct  approach  is  to  simply  retain  the  quadratic  terms  in  the 
expansion  (2.43)  and  to  express  the  scattered  fields  in  terms  of 
Fresnel  integrals.  This  latter  approach  has  in  fact  been  pursued 
at  TSC .  The  results  of  our  work  tc  date  are  presented  in 
Appendix  A. 

In  spite  of  these  various  restrictions  upon  the  applicability 
of  Equations  (2.47)  and  (2.48),  they  cam  provide  a  great  deal  of 
useful  information  about  scattering  phenomena.  In  the  next  sec¬ 
tion,  these  equations  will  be  used  to  treat  the  problem  of  locali¬ 
zer  signal  scattering  by  a  flat,  vertical  wall. 

2.5  LOCALIZER  SIGNAL  SCATTERING 

To  illustrate  the  techniques  developed  in  the  preceding 
section,  we  will  treat  the  important  problem  of  localizer  signal 
scattering  by  a  flat,  vertical,  rectangular  wall.  The  geometry  of 
the  problem  is  depicted  in  Figure  2.5.  The  x-y  plane  is  the 
ground  plane.  For  simplicity,  it  is  assumed  that  the  ground  is 
perfectly  flat  and  perfectly  conducting.  The  x-axis  represents 
the  center  line  of  the  runway.  The  z-axis  is  perpendicular  to  the 
ground  and  points  out  of  the  page.  The  unit  vectors  in  the  x,  y, 
and  z  directions  are  denoted  by  ex,  e^,  and  ez  respectively.  The 
localizer  antenna  is  located  a  distance  H  above  the  ground  and  has 
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coordinates  (0,0,H)  relative  to  the  x,  y,  z  coordinate  system.  H 
is  typically  on  the  order  of  10  feet  for  the  localizer. 


Figure  2.5.  Scattering  From  a  Vertical  Wall, 

Ground  Plane  Shown 

In  the  far  field  of  the  localizer  antenna,  the  electric 
field  can  be  written  as  follows: 

x)cR 

E  =  G  Eq  f($)  ~ —  ,  (2.50) 

where  u  is  a  unit  polarization  vector,  EQ  is  an  amplitude,  R  is 
the  distance  from  the  localizer  to  the  field  point,  and  f(v)  is 
the  horizontal  field  pattern  of  the  antenna.  For  an  arbitrary 
field  point  with  coordinates  (x,y,z),  R  and  $  are  defined  as 
follows: 

^x2+y2+(z-H)2 

. _  “1  (y  \  II  Cll 


R  = 


. . . . . . Hi . . . . . 


The  electric  field  generated  by  the  localizer  is  assumed  to  be 
polarized  parallel  to  the  ground  so  that  u  is  given  by 


ye  -  x  e. 
x  x  j 


(2.52) 


VHS — J 
x  +y  i 
p 

the  field  point. 


+y  is  the  horizontal  distance  from  the  antenna  to 


To  account  for  reflections  from  the  perfectly  conducting 
ground  plane,  an  image  antenna  located  at  (0,0,-H)  is  included 
in  the  formulation  to  cancel  the  E  field  of  the  antenna  on  the 
ground  (z=0) .  The  total  incident  electric  field  (direct  plus 
ground  reflected)  is  given  by 


E.  =  u  E  f($) 
r  o 


(2.53) 


Vr  o  2  2 

x  +y  +{z+H)  is  the  distance  from  the  image  antenna  to 
the  field  point. 

Equation  (2.53)  can  be  rewritten  as  follows: 


Ef  =  u  E0  f(«) 


ikR  f  Reik{Ri-R)l 

V-  1  -R7 - —  I*  (2*! 


We  now  make  the  assumption  of  small  elevation  angles.  That  is  to 

Vr_2  2 

x  +y  is  much  greater  than  the 

magnitudes  of  (z-H)  and  (z+H) .  If  Dp>>|z-H|  and  |z+H|,  then 
Dp=R~R^  and  the  factor  R/R^  in  Equation  (2.54)  can  be  set  approxi¬ 
mately  equal  to  1.  However,  the  path  length  difference  (R^-R)  in 
the  exponential  must  be  approximated  more  accurately  because  of 
the  oscillatory  behavior  of  the  exponential.  To  this  end,  we 
expand  R  and  R^  in  power  series  retaining  only  the  two  leading 
terms: 


R  =  D  + 
P 


(z-H) ‘ 
2°P 


(2. 55) 


*i  =  Dp  + 


(2.56) 


Consequently,  the  path  length  difference  (R^-R)  is  given  approxi¬ 
mately  by 


r  -  R  *  22H 

Ri  Dp 


(2.57) 


Incorporating  these  various  approximations  into  Equation  (2.54), 
we  obtain  the  following  expression  for  at  field  points  with 
small  elevation  angles: 


Si  =  “  E0 


JkR  r 


2ikzB/D 
e  P 


(2.58) 


The  corresponding  incident  magnetic  field  can  be  derived  from 
Equation  (2.58)  by  recalling  that  E  and  H  have  the  following 
asymptotic  relationship  in  the  far  field  of  a  localized  source 
{Eq.  (2.19)): 


■  w 


1/ 2  -  - 
(rxE) 


(2.59) 


Consequently,  is  given  by 


le\1/2  eikR  r 
Hi  =  (^)  Eo  f(*>  V-  - 


2ikzH/D 

e  F  (Rxu)  ,  (2.60) 


where  R=R/R  and  R=xe+ye  +(z-H)e  .  The  magnetic  polarization 

x  y  z 

vector  (Rxu)  is  given  by 


_  [yy»J 


RXU 


D  ( xe  +ye„ 

mi 

p 


(2.61) 


The  magnitude  of  the  second  term  in  Equation  (2.61)  is 

(z-Hi/R  =  |  z— H  j  /D^  which  we  are  assuming  is  very  small  compared  to 

1.  The  magnitude* of  the  first  terra  on  the  other  hand  is 
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approximately  equal  to  1  when  D^>> | z-H| .  Consequently  R  x  u 
approximately  equals  (-ez)  for  small  angles  of  elevation.  Making 
this  approximation.  Equation  {2.60)  can  be  rewritten  as  follows: 


H-  =  -  e 
x  z 


E0  f(0) 


1  - 


2ikzH/D 


(2.62) 


It  will  be  assumed  that  for  any  point  on  the  vertical  wall, 

D^>> j z— K |  so  that  Equation  (2.62)  can  be  used  to  accurately 
describe  the  variation  of  iL  over  the  surface  of  the  scatterer. 

The  length  and  height  of  the  vertical  wall  depicted  in  Figure 
2.5  will  be  denoted  by  L  and  h  respectively.  The  coordinates  of 
the  midpoint  of  the  base  of  the  wall  will  be  deonoted  by  (x^ ,  y^, 
o) .  If 

v  ■ 

is  much  greater  than  L,  then  the  angle  subtended  by  the  wall  at 
the  origin  of  coordinates  will  be  very  small.  Assuming  that  this 
is  the  case,  f(d)  in  Equation  (2.62)  can  be  replaced  by  f  (p)  where 
0=tan~2 (y^/x^)  to  a  very  good  approximation.  In  other  words,  we 
can  ignore  the  variation  of  the  horizontal  antenna  pattern  over 
the  surface  of  the  wall  when  Dd^>>L,  assuming  of  course  the  f{$) 
does  not  vary  too  rapidly  with  i- . 

Let  denote  the  distance  from  the  antenna  to  the  mid  point 
of  the  base  of  the  wall: 


R1  = 


R  +  y?  +  h2 


(2.63) 


When  R,  is  much  greater  than  the  dimensions  of  the  wall,  the  dis- 
tance  R  from  the  antenna  to  any  pci  .t  P  on  the  wall  can  be  repre¬ 
sented  approximately  as  follows: 


R  —  R^  +  •  r 


(2.64) 


where  Rj=R^/R^,  Rl=xl®x+^l®y~H®z '  311(2  r  a  vect°r  in  the  plane  of 
the  wall  drawn  from  the  mid  point  of  the  base  to  the  point  P. 
Equation  (2.64)  results  from  expanding  R=jR^+r^  to  terms  of  first 
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order  in  r.  Let  rj  denote  a  unit  vector  parallel  to  the  wall  and 
parallel  to  the  x-y  plane  (Figure  2.5).  Using  the  point  (x^,y1#0) 
as  a  new  origin  of  coordinates,  the  vector  r  can  be  represented  as 
follows: 

r  =  n  n  +  2  e2  ,  (2.65) 


where  r,  is  variable  ranging  from  -L/2  to  +L/2  and  2  is  just  the 

elevation  of  the  point  P  on  the  wall  above  the  ground.  The-  vector 

n  can  be  expressed  in  terms  of  e  and  e  as  follows: 

x  y 

n  =  cos  5  e  -  sin  S  e  (2.66) 

*  y 

where  5  is  the  angle  between  n  and  the  x-axis  (Figure  2.5).  Sub¬ 
stituting  (2.66)  in  (2.65),  we  obtain  the  following  expression 
for  r: 


r  =  n  cos  3  e^.  -  n  sxn 


+  2  e 


(2.67) 


The  path  length  difference  ft-^-r  cam  now  be  calculated: 


(2.68) 


Since  R^D^  (H«Dd1) ,  x1/R1=x]l/Dd1=  cos  i>  and  y1/R1=y1/Dpl=  sin  0. 
Making  these  approximations,  Equation  (2.68)  becomes: 

R,  •  r  =  n cos (8+0)  -  f-  .  (2.69) 

Pi 

Substituting  Equation  (2.69)  into  Equation  (2.64),  we  obtain  the 
following  approximate  expression  for  the  distance  R: 

zH 

R  =  R  +  r,  cos  (5+0)  -  jr —  .  (2.70) 

1  Dpl 

In  the  denominator  of  Equation  (2.62)  R  can  be  simply  replaced 
by  R^  and  in  the  factor  [1  -  e2i‘<zH^Dp]  ,  can  be  replaced  by  Dp^. 
Incorporating  all  of  these  approximations  into  Equation  (2.62),  we 
obtain  the  following  expression  for  on  the  surface  of  the  wall: 
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H.  = 

1 


1/2 


ikR^ 


Eo  fit) 


ikncos  (6+<0 


2i 


2i  sin  k 


(2.71) 


The  expression  for  given  in  Equation  (2.71)  can  nt,rf  be 
substituted  into  Equation  (2.48)  to  calculate  the  scattered  elec¬ 
tric  field  Eg  produced  by  the  wall.  The  origin  of  coordinates  for 
this  integration  will  be  the  aid  point  of  the  base  of  the  wall 

(xl,yl'°) * 

Equation  (2.48)  is  reproduced  here  for  ease  of  reference 


V*’) 


-ik 


(r*  -r)  ds. 


(2.48) 


The  unit  normal  vector  n  points  into  S+ (Figure  2.5).  In  terms  of 

the  unit  vectors  e  ,  e  ,  and  e  ,  n  is  given  by 

x  y  z 

n  =  sin  8  e  +  cos  8  e  (2.72) 

*  y 

Clearly,  n  it  is  parallel  to  n  since  is  parallel  to  ez  (Eq. 
(2.71)). 

The  vector  r*  is  drawn  from  the  point  (x.,y^,o)  co  the  ob¬ 
servation  point  which  is  of  course  the  location  of  the  receiving 
antenna  on  the  aircraft-  Let  x^ ,  yj,  z^  denote  the  coordinates 
of  the  receiver  relative  to  the  x,  y,  and  z  axes  respectively. 
Then  r *  is  given  by 


r’  =  (x2-xx)  ex  +  (y2~y1>  %  +  z2  *z  * 
The  magnitude  r*  of  the  vector  r*  is  denoted  by  R^: 


=  «/(x_-x.)2  + 


R2  =  V 


2  1' 


(y2-yi>2  +  z2 


(2.73) 


(2.74) 


The  horizontal  distance  of  the  receiver  from  the  point  (x^,yj,o) 


is  denoted  by  D 


p2‘ 


Dp2  =  V,X2‘X1’2  *  (y2'yli 


(2.75) 
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It  is  assuned  that  the  elevation  angle  of  r*  is  very  snail.  That 
is  to  say,  it  is  assuned  that  D^2>>z2‘  T*e  piiase  tern  r’*r  appear¬ 
ing  in  the  integral  in  Equation  (2.43)  is  given  by 


£ ' -r  =  n  cos 


( x2~x x) 


-  n  sxn  3 


ly2~yl] 


(2.76) 


where  use  has  been  nade  of  Equation  (2.67).  Since  D  _>>z_,  D  _=»  . 

pz  Z  p2  2 

Consequently,  the  following  approx inate  relationships  hold: 


(Xj-Xj)  (X^-Xj^) 

~%r 


cos  Y 


{y2_yl}  iy2~yl) 


sin  y 


(2.77) 


where  v  is  the  angle  between  the  projection  of  the  vector  r  *  onto 
the  x-y  plane  and  the  x-axis  (Fig.  2.5).  Consequently,  Equation 
(2.76)  can  be  rewritten  as  follows: 


r'-r  =  n  cos  (y~- )  +  — ~  *  (2-78) 

P2 

The  integral  in  Equation  (2.48)  can  now  be  evaluated  quite 
easily.  It  should  be  noted  that  the  range  of  integration  of  the 
variable  z  is  -  h  <  z  <  h.  This  range  of  integration  mist  be 
used  in  order  to  take  into  account  the  inage  of  the  wall  in  the 
ground  plane.  Carrying  out  the  integrations  in  n  (-L/2<ri<+L/2) 
and  z  (~h<z*+h) ,  we  obtain  the  following  expression  for  the  total 
scattered  field  Eg  (direct  plus  ground  reflected) : 


where  sinc(x)  =  sin  x/x. 


For  the  low  angles  of  elevation  which  we  are  assuming,  the 
triple  vector  product  r'  x  (r'xn)  can  be  simplified  greatly.  We 
first  expand  the  triple  vector  product  as  follows: 


r'  x  (r’xn)  =  (r'*n)  r!  -  n 
The  unit  vector  r'  is  given  by 

z. 


(x,-x- ) 

•  -  — t — i —  e 


r  = 


+  (y2~?V  g  +^e 

h2  y  r2  2 


(2.80) 


(2.81) 


For  small  angles  of  elevation  (R2=Dp2>>z2^ •  we  can  ignore  the  z- 
component  of  r’  and  use  Equation  (2.77)  to  obtain  the  following 
expression  for  r'. 


r'  =  cos  y  e 


-sin  y  e 


(2.82) 


Consequently,  r'*n  3  cosG  cosy  +  sinQ  siny  =  cos (y-0)  (Eq.  (2.66)). 
Making  these  approximations  in  Equation  (2,20),  we  obtain  the 
following  approximate  expression  for  the  polarization  vector 
r '  x  (r'xn) : 


r '  x  (r'xn)  =  -  sin 


n  (y-6)|si 


sin  ye  +  cos  y  e 

*  y 


] u  ■ 


83) 


which  is  of  course  parallel  to  the  ground  plane.  The  magnitude  of 
the  vector  given  in  Equation  (2.83)  is  just  sin  (y-0)  .  Conse¬ 
quently,  the  magnitude  of  the  scattered  field  is  given  by 


Es  =  aTTL<VW),£ 


iMRj+Rj) 


R1R2 


sin  (y-0)  sine  $=■  cos(0+i|<)  -  cos  (y-0 


'J 


sine 


kh fc-y  ' sinc  [kh [4 


(2.84) 


For  nn  omnidirectional  receiving  antenna,  the  expression  given  in 
Equation  (2.84)  is  the  appropriate  one  to  use  for  calculating 
perturbations  in  the  difference  in  depth  of  modulation. 
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2.6  COMPARISON  WITH  PREVIOUS  WORK 

The  result  obtained  in  the  preceding  section  for  the  scattered 

electric  field  produced  by  a  flat  vertical  wall  will  now  be  compared 

with  the  solutions  to  the  same  problem  obtained  by  the  Ohio  Univer- 

4 

sity  and  I.B.M  researchers. 

Our  Equation  (2.84)  for  Eg  does  not  differ  significantly  from 
the  expression  obtained  by  the  Ohio  University  group.  This  comes 
as  no  surprise  since  our  methods  of  approach  are  almost  identical. 
Their  solution  is  somewhat  simpler  in  that  it  does  not  contain  the 
sine  functions  which  depend  upon  the  height  of  the  localizer  antenna 
(H) ,  the  height  of  wall  (h) ,  and  the  height  of  the  receiver  (z^) 
that  appear  in  Equation  (2.84).  Instead,  their  expression  for  E 

3  s 

is  simply  proportional  to  h  H  z^.  The  fact  that  the  Ohio  University 
solution  has  a  simpler  dependence  upon  the  parameters  h,  H,  and  z^ 
than  does  our  solution  is  due  to  the  fact  that  the  Ohio  University 
researchers  used  a  linear  approximation  to  describe  ground  reflec¬ 
tions.  Specifically,  instead  of  multiplying  the  antenna  field  pattern 

2  x  k  z  /D 

by  the  factor  [1  -  e  '  p]  to  account  for  ground  reflections  (Eq. 

(2.58)),  as  we  did,  the  Ohio  University  researchers  multiplied  the 

antenna  field  pattern  by  (2kzH/D  ) .  That  is,  they  assumed  that 

P  2ikzH/D 

(2KzH/Dp)  <<  1,  and  essentially  expanded  [1  -  e  '  p]  in  a  power 

series  in  (2kzH/D  )  retaining  only  the  first  two  terms.  They  made 
P 

a  similar  approximation  to  describe  the  reflections  of  the  scatter¬ 
ed  fields  from  the  ground  plane.  Our  expression  for  Eg  in  fact  re¬ 
duces  to  the  Ohio  University  expression  if  the  sine  functions  in¬ 
volving  H,  h,  and  z2  appearing  in  Equation  (2.84)  are  approximated 
by  the  first  two  terms  in  their  power  series  expansions.  In 
essence  then,  there  are  no  fundamental  differences  between  our 
treatment  of  the  scattering  problem  and  the  Ohio  University  treat 
ment . 

On  the  other  hand,  the  I.B.M.  approach  to  the  problem  of 
localizer  signal  scattering  by  a  vertica)  wall  differs  fundamentally 
from  our  approach  and  that  of  Ohio  University.  For  convenience, 
their  approach  to  the  problem  will  be  illustrated  using  our  equa¬ 
tions  and  notation. 

For  ease  of  reference.  Equations  (2.28)  and  (2.29)  are 
reproduced  here: 
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(2.28) 


E  (r 
s 


i(uu  (nxH  )4*  +  (nxE  Jx^+ta-E 
s  s  s 


ds. 


Hs  (r 


■'-hf  [i 


iwe (nxE  )?-(nxH  (n-H  }$V 

3  S  S 


ds 


(2.29) 


Equations  (2.28)  and  (2.29)  express  the  scattered  fields  outside 
a  scatterer  in  terms  of  integrals  of  the  scattered  fields  over  the 
surface  of  the  scatterer.  For  treating  perfectly  conducting 
scatterers,  we  chose  to  work  with  Equation  (2.29)  since  the  inte¬ 
gral  simplies  greatly  for  a  perfect  conductor.  The  integral  in 

(2.29)  was  first  modified  by  adding  to  the  right  hand  side  of 

(2.29)  the  corresponding  integral  of  the  incj.dent  fields  E.  and  It. 
Since  the  surface  integral  in  (2.28)  and  (2.29)  of  it  and  It  over 
S1  are  identically  zero,  the  left  hand  side  of  (2.29)  was  un¬ 
changed  by  this  addition  and  we  obtained  the  following  modified 
expression  for  Hs(r'): 


Hg(r')  =  I7  J  iue(nxE)y-(£xH)x7?-(n-&$¥  ds,  (2.31) 


where  E  and  H  are  the  total  fields  on  S. (E=E.+E  ,  H=H.+H  ) . 

1  x  s  is 

The  analogous  expression  for  the  scattered  electric  field  is: 


■  ■  W  [* 


iop  (nxH)'P  +  (nxE)  x$¥+ 


(n-E)^j 


ds 


(2.85) 


For  a  perfectly  conducting  scatterer,  the  surface  is 
divided  into  an  illuminated  side  S+  and  a  shadow  side  S_.  On  the 
shadow  side,  it  is  assumed  as  a  first  approximation  that  2h5=0 
(Eg=-E^  and  Hs=-H^  on  S_) .  With  this  approximation,  the  integrals 
in  (2.31)  and  (2.85)  are  extended  only  over  the  illuminated  side 
S+  of  S-^.  That  is  to  say,  Es(r’)  and  Hs(r*)  are  expressed  in 
terms  of  surface  integrals  of  the  total  fields  over  S  . 

*T* 
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Had  we  chosen  to  work  directly  ich  Equations  (2.28)  and  (2.29) 
rather  than  (2.31)  and  (2.85)  the  same  result  would  have  been  ob¬ 
tained.  For  example,  assuming  th  Eg  =  -  E^  and  Hg  =  -  on  S_, 
Eouation  (2.28)  becomes: 


Es(r’> 


‘ '  hi  [iu“ 


(nxH  )y+(nxE  )x7f+ 

5  S 


f+(n*Es)V?J< 


J  J^iwu  (nxH^)'r  +  {nx£^)x^?  +  (n-iL  )?yj( 


(2.86) 


However,  we  know  that  the  surface  integral  of  E^  and  over  the 
total  surface  is  identically  zero.  Consequently,  the  surface 
integral  in  (2.86)  of  E.  and  5.  over  S_  is  the  negative  of  the 
corresponding  surface  integral  of  E^  and  over  S+.  Therefore, 
Equation  (2.86)  can  be  rewritten  as  follows: 

Eg (r*)  ~  ~  J  iwu(nxHs)y+(nxis)x^¥+(n*Es)^¥  ds 
'  s+  L 

— J  J^it (nxH^)'?  + (nxE^)x7?+ (n*E^)  ds 

S+  (2.87) 

Combining  these  two  integrals  over  S  and  recalling  that  E=E.+S 
■+  +  1  s 
we  are  left  with  a  surface  integral  over  5+  of  the  total 

fields  E  and  H. 

The  I.B.M.  researchers  essentially  use  the  same  integral  Equa¬ 
tions  (2.28)  and  (2.29)  to  treat  the  problem  of  localizer  signal 
scattering.  To  obtain  an  expression  for  the  scattered  electric 
field  produced  by  a  vertical  wall,  they  simply  carry  out  the 
integral  appearing  in  (2.28)  over  the  illuminated  surface  cf  the 
wall.  That  is  to  say,  they  use  the  following  expression  for  Es(r‘): 
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VJ,>  -  -  hf  [iul 


(nxHs )  ?  +  (nxEg )  x%+  (n  •  Eg ) 


ds  . (2.68) 


That  is  to  say,  instead  of  integrating  the  total  fields  over  S+ 

as  we  did,  they  just  integrate  the  scattered  fields.  In  effect. 

they  ignore  the  integral  of  E  and  H_  over  S_,  the  shadow  side  of 

S^.  Instead  of  setting  £s  and  Hg  equal  to  — E^  and  -H^  respectively 

on  the  shadow  side  S  ,  they  in  effect  assume  that  E  and  H  are 

—  s  s 

identically  zero  on  S  which  is  not  true. 


Comparing  our  result  for  the  magnitude  Es  of  the  scattered 
field  produced  by  a  vertical  wall  with  the  I.B.M.  results,  we  find 
that  the  only  real  difference  between  the  two  expressions  is  that 
whereas  our  Equation  (2.34)  contains  a  factor  [2  sin  { y— 6 ) 3  co 
describe  the  horizontal  field  pattern  cf  the  wall,  their  Equation 
(2.43)  contains  a  corresponding  factor  MG  which  is  given  in  our 
angle  notation  by: 


■[- 


sin(9+$)+sin(Y 


-9)  cos  y-'i—2v  I  . 

J  L  •• 


(2.89) 


The  effects  of  this  difference  on  perturbations  in  the  difference 
in  depth  of  modulation  predicted  by  the  two  models  can  be  very 
pronounced  as  can  be  seen  in  some  of  the  accompanying  comparative 
graphs.  Figures  2.6  and  2.7  (static  and  dynamic)  show  the  dif¬ 
ference  between  the  new  and  old  (I.B.M.)  formulations.  Also,  as 
examples  of  the  dependence  of  the  derogation  on  the  building  size, 
location  and  orientation,  we  show  in  Figures  2.8  to  2.18  the  micro¬ 
amp  deviation  of  a  localizer  V-ring  signal  that  would  be  received 
by  an  aircraft  flying  the  centerline  of  the  runway  in  level  flight, 
50  feet  above  the  ground  as  predicted  by  the  nev?  math  model  for 
scattering  off  a  vertical  wall. 

These  Cal  Comp  generated  figures  are  presented  hero  as  ex¬ 
amples  of  typical  output  obtained  from  the  computer  model  and  as 
examples  of  one  way  in  which  this  output  cai\  be  presented.  In  all 
figures  the  x-distance  is  measured  from  the  Localizer <as  origin 
and  the  y-distance  from  the  centerline  of  the  runway. 
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In  Figure  2.8,  the  microampere  deviation  calculated  along  the 
runway  centerline  caused  by  a  100  foot  by  50  foot  rectangular  wall 
situated  behind  the  Localizer  and  470  feet  from  the  runway  center- 
line  is  shown  to  increase  markedly  as  the  wall  is  moved  closer  to 
the  Localizer.  This  is  expected  as  more  energy  is  reflected  back 
onto  the  runway  in  these  cases. 

Figure  2.9  shows  a  similar  marked  increase  in  the  derogation 
for  buildings  located  closer  to,  but  in  front  of  the  Localizer. 

Figures  2.10  and  2.11  show  the  derogation  due  to  a  larger 
building,  500  feet  by  50  feet,  which  is  located  further  from  the 
runway  centerline  (1000  feet  instead  of  470  feet)  .  Well  beyond 
the  Localizer  in  either  direction,  (-2000  feet  or  +8250  feet), 
there  is  only  negligible  derogation,  while  closer  in,  at  -1000 
feet  and  at  +500  fee:,  the  derogation  is  significant. 

The  explicit  variation  of  the  derogation  on  the  building's 
position  relative  tc  tne  runway  centerline  is  shown  in  Figures 
2.12  and  2.13.  There  is  a  dramatic  increase  in  the  amount  of 
derogation  when  the  building  is  placed  close  to  the  runway  center¬ 
line  (e-g.,  when  y=4  J  feet). 

Figures  2.14  and  2.15  show  the  derogation  for  different 
orientations  of  the  building,  bctn  in  terms  of  the  magnitude  and 
the  location  of  the  derogation  on  the  runway  centerline.  As  in¬ 
dicated  in  the  accompanying  sketches,  the  different  locations  of 
the  derogation  along  the  runway  centerline  occur  because  of  the 
different  specular  reflection  directions  of  the  various  building 
orientations. 

Finally,  in  Figures  2.16-2.18,  the  increase  in  the  derogation 
with  the  increase  of  building  size  is  shown. 
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angle  alpha  is  varied  to  givo  specu 
parallel  to  the  runway. 
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atant.  The  angle  alpha  is  varied  to  give  specular 
rofloction  parallel  to  the  runway. 
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200ft.  are  used  while  the  "X"  and  "5f"  positions 
and  the  height  of  the  wall  are  held  constant. 


3.0  MODELING  OF  ILS  SIGNAL  DETECTION 


3 . 1  INTRODUCTION 

In  implementing  computer  prediction  of  the  course  deviation 
indication  (CDI)  for  complex  derogating  environments  and  arbitrary 
selection  of  the  localizer  and/or  glide  slope  signal  transmitting 
system,  it  is  desirable  to  have  a  unified  approach  to  the  modeling 
of  DDM  detection  by  standard  aircraft  ILS  receivers-  To  effec¬ 
tively  treat  any  existing  or  future  ILS  system,  a  detection  model 
must  adequately  account  for  systems  utilizing  one  or  two  carrier 
frequencies,  for  arbitrary  relative  phasing  between  different 
signal  components,  variation  of  receiving  antenna  gain  patterns 
and  effects  of  aircraft  speed.  The  models  presently  in  use  fall 
short  of  this  objective.  The  X.B.M.  code  is  developed  only  for 
the  localizer  null  reference  system;  the  Ohio  glide  slope  programs 
incorporate  null  reference,  sideband  reference  and  capture  effect, 

5 

but  in  the  latter  case  the  auxrlliary  carrier  signal  is  ignored. 

The  Ohio  codes  are  further  restricted  in  applicability  to  moderate 
derogations  near  or  on  the  glide  course  centerline,  and  none  of 
the  foregoing  models  take  any  account  of  relative  Doppler  effects. 

Below  we  develop  a  unified  model  of  ILS  signal  reception  that 
attempts  to  reach  a  higher  level  of  completeness  and  correctness 
for  actual  receivers  than  the  previously  used  models.  In  the 
course  of  the  derivation,  a  number  of  approximations  are  inevitably 
taken.  While  it  is  felt  that  these  are  all  reasonably  justified, 
the  overall  consequences  of  the  model  should  be  verified  by  empir¬ 
ical  data.  The  treatment  does  illuminate  the  relative  signifi¬ 
cance  for  the  detected  ILS  signal  of  carrier  and  sideband  dephas¬ 
ing  and  of  Doppler  effects  on  receiver  dynamic  response. 

3.2  RECEIVER  MODEL 

Localizer  and  glide  slope  signals  are  assumed  to  be  detected 
by  heterodyne  AM  receivers  similar  in  basic  design  to  the  sche¬ 
matic  diagram  shown  in  Figure  3.1.  The  figure  indicates  generally 
how  the  90  to  150  Hz  modulations  are  amplified,  detected,  and 
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separated  and  the  difference  applied  to  one  of  the  movements  of  a 
cross  pointer  indicator.  The  sum  of  the  90  and  150  Hz  outputs  is 
held  constant  by  an  AGC  feedback  loop,  so  that  once  calibrated,  the 
instrument  gives  a  faithful  readout  of  tne  apparent  CDI  {course 
deviation  indication)  whenever  the  incident  signal  radiation  ex¬ 
ceeds  a  minimum  level. 

The  guidance  radiation  field  (localizer  or  glide  slope)  is 
coupled  to  the  receiver  input  via  an  antenna  which  may  be  assumed 
to  have  specified  directional  and  polarization  characteristics. 

Any  antenna  used  in  practice  may  be  treated  as  a  linear  system. 

This  fact  allows  us  to  consider  the  scalar  current  transmitted  to 
the  receiver  input  load  as  the  resultant  of  individual  currents 
generated  by  various  electromagnetic  field  components  incident  on 
the  antenna.  We  shall  assume  that  all  incident  component  fields 
are  transverse  and  harmonic,  and  so  are  adequately  described  by  an 

A 

electric  vector  E  and  a  unit  vector  k  in  the  direction  of  propaga¬ 
tion.  For  a  number  of  such  component  fields  incident  on  a  speci¬ 
fied  antenna,  the  receiver  input  current  can  be  expressed  in  the 
form: 


=  Z  gtkj)  •  E.  (3.1) 

where  g (k)  conveys  both  the  gain  and  polarization  characteristics 
of  the  antenna  for  plane  waves  incident  in  the  direction  k. 

The  ILS  radiation  fields  to  be  considered  may  be  generated  by 
multiple  element  antennae  and/or  waveguide  radiators  operating  on 
one  or  two  carrier  frequencies.  In  every  case,  a  number  of  dis¬ 
tinct  directional  modulations  are  radiated,  each  in  a  specific 
spatial  pattern.  The  separate  directional  waveforms  are  comprised 
of  certain  combinations  of  90  Hz  and  150  Hz  amplitude  modulations 
with  prescribed  carrier  phases,  added  to  unmodulated  carrier  in 
fixed  ratios.  The  measured  CDI  at  a  particular  aircraft  location 
is  based  on  the  relative  strengths  of  90  Hz  and  150  Hz  signals 
detected  in  the  resultant  of  all  direct  and  scattered  ILS  radia¬ 
tion  reaching  the  receiving  antenna. 
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tern  Fourier  series: 


S„<t>  =  S=3  K3 (t)  -  S„5  U5(t) 


sA(t>  =  s-  u_(t)  +  s*  u.(t). 


(3.3) 


ire 


U_  (t)  =  cos  60 


The  terss  C^,  ,  S  S.'_ 


(3.4) 

include  as  factors  the  relative 
strengths  of  the  initially  radiated  carrier  and  sideband  amplitudes, 
the  path  attenuation  due  to  field  divergence,  reflections,  etc., 
and  the  earn  factor  for  the  direction  of  incidence  on  the  receiving 
antenna  g(k  } *E  .  The  total  complex  IF  output  signal  is  the  sum 

D  O 


Souaticn  (3.2)  over  all 


It loath 


CCS  o . 


This  is  passed  to  the  second  detector  which  generates  an  audio 
frequency  signal 

VAptt)  =  i  I  VIFo(t)|.  (3.5) 

P 

Determination  of  the  CDI  is  based  on  the  relative  amplitudes  of 
90  Hz  and  150  Hz  components  found  in  this  audio  signal  by  passing 
it  through  selective  filters.  Denoting  the  action  of  these 
filters  symbolically  by  the  operators  Hg0  and  (to  be  defined 

presently) ,  we  may  represent  the  outputs  of  an  ILS  receiver  as 


V90(t' 

-  "so  lvir(t>1 

(3.6a) 

V150(t) 

=  H150  tyAr{t,1i 

(3.6bi 

CDI  =  S 

V  -V 

150  90  microamperes 

(3.7) 

m 

V150+V90 

In  Ecuation  (3.7)  the  factor  s_  is  t _ .  cross-pointer  sensitivity 

which  has  the  standard  values  387. C  for  localizer  (horizontal) 
signals  and  576.0  for  glide  slope  (vertical)  signals.  The  magni¬ 
tudes  V, ana  V_n  are  obtained  bv  rectification  of  the  detected 
laO  90 

modulation  tone  signals  and  (t) . 

3.3  AUDIO  SIGNAL  DETECTION 


3.3.1  Capture  Effect 


we  now  seek  a  reasonably  accurate  estimate  of  the  detected 
90  Hz  modulations  in  the  general  multipath,  capture  effect  situa¬ 
tion.  for  this  we  need  an  expression  for  V._{t)  in  the  form  of  a 

*  r%£ 

Fourier  series,  and  we  need  suitable  definitions  of  the  filter 
operators  H9Q  and  H150- 


It  is  the  possible  presence  of  many  signal  derogating  scatter 
ers  in  the  glidepath  environment  that  makes  analytical  derivation 
of  V._(t)  from  V _ ,(».'•  difficult.  Because  of  the  receiving  air- 

nr  i.r 

craft's  motion,  each  derogating  signal  is  received  with  a  possibly 


significant  Doppler  shift  : 


ILS  transmission . 


relative 


to  the 


The  consequence  is  that  a  range  of  frequencies  may  generally  be 
present  in  the  audio  output  of  the  second  detector,  instead  of  sole¬ 
ly  the  modulation  tone  frequencies.  The  action  of  the  90  and  150 
Hz  filters  on  such  complex  signals  must  be  considered. 

Additional  spurious  signals  occur  in  reception  of  radiation 
from  capture  effect  transmitting  systems  near  all  harmonics  of  the 
8  kHz  carrier  separation  frequency.  This  fact  allows  a  convenient 
reduction  of  the  audio  signal  analysis  problem,  if  the  capture 
effect  of  the  detector  is  treated  first,  lie  follow  a  procedure 
suggested  in  part  by  Manney.®  The  total  IF  output  signal  may  be 
written  as: 


Vjplt)  = 


=  i 


A1(t) 


i“ct 
e  c  A- 


(t) 


— lta_t 
C 


15.8) 


where: 


A1(t)  =  X  (C  +  So)  e1^ 
P 

V 

A2 (t )  =  L  (c;  +  S')  e  p 

n  *  * 


a  =  a  -  a' 
C  c  c 


(3.9) 


The  audio  signal  (Eq.  (3.5))  can  now  be  expressed  as 


VAF(t) 


=  h 


A1(t)i2  + 


|A2(t) 


-r  2  Re  ^ 


(t)A_ (t)*e 


-it  t 


2 

plex  frequencies  A. 


)]•' 


10) 


Now  the  scale  of  the  time  variation  of  the  cc 
and  A^  is  given  by  the  modulation  frequencies  90  Ez  and  150  Hz  or 
perhaps  by  the  latter  frequency  shifted  by  the  upper  limit  of  ob¬ 
served  Doppler  displacement.  In  the  worst  case  for  glide  slope 
signal  reception,  the  frequency  spectral  range  of  A.  and  A^  is 
certainly  less  than  500  Hz,  which  can  be  considered  small  compared 
to  & c/2x  =  8  kHz.  Accordingly,  we  may  expand  the  instantaneous 


value  of  V„„  (t)  in  a  cover  series  in  exo  (io  t) 

nr  w 


X.e.  : 


v!i  =  2  > 


nii. 


VR(t)e 


(3.11) 
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In  this  expression,  because  of  the  spectral  narrowness  of  the 
factors  ( t ) ,  the  term  primarily  contributing  frequencies  within 

the  passbands  of  the  modulation  filters  is  VQ(t).  This  term  can 
be  defined  by  averaging  the  expression  for  V„„(t)  (Eq.  (3.10)) 
over  the  phase  of  exp  (iu^t) ,  i.e.t 

2- 


V0(t)  =  1/2 


/ 


Ax(t)i2+  |A2(t)  i2  +2Re(A1A 


*e-)J 


11/2 


dp  (3.12) 


This  integral  is  easily  transformed  to: 

s/2 


v0(t)  = 


where : 


2(iA1(t)l  4-  |A2(t)})/-  £ 


2  2  ^2 
(1-k  sin  5)  ds 


k2  =  4jA1(t)A2(t){/(|Ai(t)  j  -r  !A2(t)i)2 


a  =  (o+512)/2 


io. 


(3.13) 


e  12  =  Aj^Aj/  ]  AjA2  i 

The  integral  on  the  right  of  Equation  (3.13)  is  a  complete  elliptic 
integral  of  the  second  kind  and  is  ccsmonly  denoted  by  £(k): 


(3.14) 


Thus  for  the  purpose  of  determining  the  detected  90  and  150  Hz 
modulations,  tlx*  unfiltered  audio  signal  is  to  a  good  approx irta tier 


V^tt)  =  |  (|A1(t)|  4-  |A.(t}|)E{k) 


(3.15) 


To  put  the  above  expression  in  usable  form,  vs  must  carry  the 
approximation  somewhat  further.  Let  the  magnitude  of  amplitude 
A^(t)  be  the  sum  of  a  constant  average  value  aQ1  and  a  zero-mean 
time  varying  function  a. (t): 

|A1(t)  |  =  aQ.  4-  a,  (t)  (3.16a) 


f 

J 
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The  quantity  a^ (t)  represents  the  total  modulation  that  would  be 
detected  if  the  second  carrier  were  completely  absent.  Similarly, 
let 


i  A2  (t)  i  =  aQ2  +  a2(t> 


(3.16b) 


Inserting  Equations  (3.16)  into  (3.15)  and  treating  a^(t)  and  a2<t) 
as  small  quantities,  we  obtain  an  expansion  for  V„ : 


lo  !o  |o  (a«)  W  (*■>)' 


V*1’  =  a0 


(3.17) 


Here: 


I  ' 


ao 

=  a01  + 

alt) 

=  (t) 

4 

=  4aoia 

_n 

j 

=  1  ;  n= 

pqr 

=  0  ;  n 

=  (-l)r 

n  >  1 

ve  find 

that  the 

+  a02  *  ' 

a][  (t)  +  , 

**n  « 2ft 

1 ooo*o 

a01+a02 

a02~a0l\ 

a02  / 

«2 

.-1)  (a01 

-a02>2/2 

2n(2n) I (r+2n-2) 1 
24n(nl)2in-p)!  p!  (n-q) 


(3.13) 


.ri!0p2\  (t) 

a01  /  * 


*01  02/ 


(a1(t)/a01-a2(t)/a02 


H 
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{3.19) 


4  3 

+  0{kJaJ) 


=  {aoi+a 

02* 

(1  a0la02/a0 

+  ax(t) 

(1 

2.  2 

-  a02  /a0  "  *** 

+  a2(t) 

2  2 

(1 

‘  aoi  /ao  -  *** 

+  ... 

Sussing  a  finite  nunber  of  the  terras  constant  and  linear  in  a^(t) 
and  a_(t)  gives  an  accurate  representation  of  the  audio  signal 
V„„(t)  when  the  parameter  kn<<l.  The  convergence  of  the  series  as 
k*l  is  not  easy  to  discuss  because  the  raodulations  of  ILS  trans¬ 
missions  are  quite  high  (40%  in  the  case  of  localizer  signals,  80% 
for  glide  slope  signals) . 


3.3.2  Decomposition  of  Single  Carrier  Envelopes 

It  reraains  to  find  suitable  Fourier  expansions  for  the  sepa¬ 
rate  carrier  envelopes  A1 (t)  and  (t ) .  In  this  section  we  will 
assurae  the  position  of  the  receiver  to  be  stationary,  so  that  the 
chases  4,5"  are  constant.  In  the  next  section  we  treat  Dcooler 

p  p 

effects  by  varying  these  phases. 


The  analysis  proceeds  below  for  A^(t)  and  is  identical 
Aj  (t) .  We  have 


where 


for 


(3.20) 

(3.21) 

(3.22) 
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We  may  put: 


S(t)e 


C  =  C  e 


r=  X  +  iY  , 


(3.23) 


(3.24) 


and,  if  |s|  <  |C|,  we  have: 


Ax(t)  =  i C | jl+2M/ i C !  +  (X2+Y2)/|C|2J 1/2 

=  !  c !  +  x  +  y2  (2  j  c  |  > _1  (i-x/  j  c  |  +  ...) 


(3.25) 


The  quantity  X  represents  the  sideband  component  of  the  received 
signal  which  is  in  phase  (or  exactly  out  of  phase)  with  the  un¬ 
modulated  carrier.  The  quantity  Y  represents  the  sideband  compon¬ 
ent  which  is  in  quadrature  with  the  carrier.  The  nonlinear  terms 
of  Equation  (3.25)  may  be  further  expanded  in  a  series  of  all  har¬ 
monics  of  tfe  fundamental  frequency  30  Hz.  However,  we  note  that 
the  lowest  order  term  Y  / (2  j  C | )  contributes  no  correction  to  the 
90  Hz  and  150  Hz  frequencies.  The  higher  order  corrections  will 
be  considered  unimportant  for  the  present  application .  Thus  in 
the  notation  used  above,  we  have 


*01  =  icl 


(3.26) 


ax(t) 


=  Z(v 

-  II  a 


3(t)  +  SppjU,.  (t)  Icos  ($^-9-^) 


p  m 


Similarly, 


a  .  U  (t)cos ($  -8. ) 
mlp  m  pl 


a02  ' C ' 


a2  (t) 


-II 

m  p 


m2p  m 


U  (t)cos ($ ’-8,) 


(3.27) 


(3.28) 


(3.29) 


In  the  above,  the  carrier  phases  are  determined  by 


=  arctan 


I  Cpn 


C  cos  (J> 

pr.  pn 


(3.30) 


where 


(Cpl,  4>pl)  correspond  to  (Cp, 

and 

(Cp2>  cp2)  correspond  to  (C^,  c^)  . 

3.3.3  Doppler  Effects 

As  a  landing  aircraft  moves  there  are  relative  changes  in  the 
phases  of  the  various  multipath  radiation  components  of  the  ILS 
signals  received.  The  amplitudes  of  the  90  Hz  and  150  Hz  terms  of 
the  detected  modulation  thus  vary  with  time,  and  accordingly,  each 
tone  acqui  s  a  frequency  spread.  Since  the  relative  Doppler  fre¬ 
quency  between  direct  and  scattered  radiation  can  be  as  much  as 
120  Hz  in  some  cases  for  typical  landing  approach  speeds,  the 
bandwidth  of  the  modulation  filters  becomes  important  in  determin¬ 
ing  receiver  dynamic  response.  The  effect  envisioned  here  is  that 
certain  greatly  Doppler  shifted  signal  components  may  either  be 
severely  damped  by  one  filter  or  unexpectedly  passed  by  the  other 
depending  on  the  circumstances  of  receiver  design,  aircraft  speed, 
etc.  This  is  independent:  of  the  intentional  RC  damping  that  is 
always  introduced  into  the  cross  pointer  indicator  circuits. 

To  evaluate  the  possible  magnitudes  of  such  phenomena,  we 
consider  a  particularly  simple  model  for  the  frequency  response 
of  each  modulation  filter.  Though  the  characteristics  of  this 
filter  may  differ  markedly  from  those  of  filters  in  actual  re¬ 
ceivers,  the  model  is  adequate  for  studying  the  general  effect  of 
filter  bandwidth  on  receiver  rejection  of  ILS  interference. 

The  two  detected  tones  which  emerge  from  the  modulation  fil¬ 
ters  were  earlier  denoted  by  V^ft)  and  v^5g  (Eq.  (3.6)).  We 
now  define  these  more  precisely:  Let 

Crt 

V90(t)  =  f  h90(T)  VAF(t-T)dT,  (3.31) 

(150)  "0  (150) 

where  h-fi)  is  the  impulse  response  function  of  a  narrowband  filter 
i  7 

corresponding  to  center  frequency  f.  Following  Middleton  ,  we 
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represent  the  narrowband  filter  by  the  response  function 

hj;(T)  =  h0(T)cos  (2iTfT+£)  (3.32) 

and  we  select  as  the  particular  "window  function" 

hQ(T)  =  2/T  0<t<T 

=  0  otherwise.  (3.32a) 

The  output  of  such  a  filter  when  driven  by  a  sinusoidal  input 
A  cos (2nf „t+o£)  is  approximately 

Vf(t)  =  (A/2)  sine  (Tr(f-f£)T)  . 

cos [2rft+6  a+ty+2v (f .-f )  (t=T/2)]  (3.33) 

JC  X* 

whenever  f+f£^>T  ^ .  The  reciprocal  of  the  intergration  time  T  is 
the  nominal  filter  bandwidth.  The  phase  is  purely  additive  and 
so  m ty  be  dropped  without  loss  of  generality.  For  brevity  we  will 
put 

Hq  (f )  =  sine  (irfT)  . 

To  derive  results  pertaining  to  the  ILS  audio  signal  (Eg.  3.5 
from  Equation  3.33,  we  assume  that  the  phase  shifts  of  individual 
multipath  components  due  to  aircraft  motion  lead  to  sinusoidal 
time  variations,  and  we  will  neglect  the  time  variation  of  ampli¬ 
tude  coefficients.  To  this  end  we  may  write  Equation  (3.19)  as 

VAF(t)  ”  B0(a01*a02>  +  Biai(t|  +  B2a2(t)  (3-34) 

where  the  B's  can  be  defined  in  terms  of  slowly  varying  amplitudes 
that  will  be  considered  constant  and  the  a's  (as  defined  by  Equa¬ 
tions  (3.26)through(3.30))  will  contain  typical  terms  such  as: 

a  U  (t)  cos  (6  -0  )  -  1/2  a  (cos  (607rmt+6  (t)-9  (t) ) 

mnp  m  np  n  mnp  np  n 

+  cos  (60-nrnt-bn  (t)  +  ©n(t)))  (3.35) 

This  form  shows  that  for  each  multipath  component  the  audio  modula 

tion  frequencies  are  modified  by  sum  and  difference  combinations 

with  the  relative  Doooler  freouenev  d/dt(o  -6  ).  Considering  all 

*  np  n 
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combinations  of  input  modulation  frequency  and  Doppler  shift  which 
may  give  net  frequencies  near  the  filter  center  frequencies  f,  we 
find  the  following  significant  terms  in  the  filter  outputs  V^(t) 
^notation  will  be  detailed  below) : 

1.  The  average  carrier  level  received  by  paths  with  lavas 
Doppler  shifts  gives  rise  to  a  contribution  from  each 
filter 


77a  H. (f+ <o  -6  J /2~) ccs (2“fT/2+S  -o  ) 
t—  nnn  0  nt  n  '  7  n  nn 


0  n  o  on?  0  '  np  "n 


(3.36) 


2.  The  modulation  at  one  frequency  modified  by  a  node  vats 
Doppler  shift  gives  rise  to  a  contribution  in  the  output 
of  the  opposite  frequency  filter 

1/2  J  BnS  a  H0(Af)cos(60irmt?(2TtAfT/2+Sn-inD))  (3.37) 

n  p  ^ 

Here  the  upper  or  lower  sign  is  taken  for  filter  frequency 
f  =  90  Hz  or  150  Hz  respectively  and 

if  =  i f-30  1  +  (0-01/2 ". 
m  np  n 

3.  In  the  case  of  equal  modulation  and  filter  frequencies 
the  output  contribution  is 


y  B  y  A  H„({s  -9  )/2r)cos2-ft 

*—  n*—  mnD  0  no  n 
n  d 


(3.38) 


The  total  detected  signal  at  each  modulation  frequency  f  is 
the  sum  of  the  contributions  1,  2,  and  3.  It  may  be  noted  that 
only  the  phase  of  3  is  well  defined  in  relation  to  the  phase  of 
the  transmitted  modulation.  This  confronts  us  with  the  final 
problem  of  this  development:  we  must  calculate  the  rectified 
signal  derived  from  each  V^(t)  which  is  passed  directly  to  the 
cockpit  panel  instruments  to  generate  the  CDI  display.  A  reason¬ 
able  approximation  can  be  reached  by  assuming  the  relative  phases 
of  contributions  1,  2,  and  3  to  be  random  and  defining  the  final 
ILS  signals  as  being  the  r.m.s.  envelope  of  the  Vf(t)  obtained 
above.  We  thus  arrive  at  the  final  formula: 


=  2  Hn(f+(.j  -a  )/2 -}■ 

0  2_,Z-  mnp  0  np  n 


m  =  0 


+  1/4  y  B2  y  a  2  •«  -  (f-30  :($-  -5  )/2a)‘ 

'  n  mnp  0  m  no  n 

n=l  d  r 


f  /  30m>0 


I,  3n  I 


n=l 


a  cos 
mnp 


Hr((i  -5  )/2r)f 
no  n  0  np  n  1 


1/2 


t  =  30mJ 
(3. 35} 


Again  the  upper  and  lower  signs  refer  to  filter  frequencies  of  90 

Hz  and  150  Hz  respectively-  As  a  practical  consideration,  each  of 

the  terms  of  Equation  (3.39)  can  be  accumulated  as  the  individual 

oath  attenuation  factors  a  are  calculated,  exceot  for  the  co- 
2 

efficients  B  ,  &  and  B  ,  the  reference  ohases  S  and  Dopoler  fre- 
.  o  n  n  n 

quencies  5^.  To  facilitate  this  method  of  computation  we  adopt 
the  approximation 


a 

wn 


dF"-' 


(3.40) 


where  R  is  the  instantaneous  range  of  the  receiving  aircraft  to 
the  transmitting  antenna  and  a  is  the  carrier  wavelength.  The 
remaining  notation  of  Equation  (3.39)  is  explained  in  the  following 
tabulation: 


f  =  modulation  filter  frequency  90  Hz  or  150  Hz. 

m  =  transmitted  harmonic  of  30  Hz;  m=3  for  9C  Hz,  n=5  for 
150  Hz,  m=0  for  constant  or  unmodulated  carrier  terms. 

0  =  Dhase  delav  of  a  oarticular  transmission  oath  to  the 

np 

receiving  antenna. 

■a  =  time  derivative  of  v 
np  np 

S  =  aooroximate  ohase  of  the  aggregate  of  direct  and  reflected 
n 

carrier  amplitudes  at  the  receiver  relative  to  the  trans¬ 
mitted  carrier  phase  (defined  by  Equation  (3.30)).  For 
course  signal  n=i  and  for  clearance  signal  n=2. 
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amno  "  carr*er  <m=0)  or  modulation  (m=3 ,5}  amplitude 
component  in  the  received  audio  signal  due  to 
transmission  by  path  p  on  the  course  (n=l)  or 
clearance  (n=2)  carrier. 


-.n 


3  =  \  r"  v 

SC  Z,  •  000  0 
n=0 


2n 


B1  2  1 000*0  -1  n(a01_a02)/a013 


B„ 


n=0 


- 1  n 


ooo*o 


2n 


^  n^a02~a01^a02^ 


n=0 


Fqoo'  =  s?ec^a*  coefficients  (defined  in  Equation  (3.18). 

The  quantities  V15Q  and  VgQ  determined  by  Equation  (3.39)  are 
inserted  in  Equation  (3-7)  to  obtain  the  undamped  CDI  which  would 
be  generated  by  a  moving  ILS  receiver.  The  “dynamic"  CDI,  which 
incorporates  the  effect  of  the  EC  damping  circuit,  may  be  derived 
from  the  undamped  CDI  by  the  usual  numerical  techniques  of  filter 
simulation. 


3.4  MODEL  STATUS 

The  detection  model  reported  here  has  been  incorporated  into 
a  computer  multipath  interference  simulation  program  which  is 
currently  under  development -  To  this  writing,  all  results  have 
been  consistent  with  expectation.  Additional  confirmation  will  be 
sought  by  comparing  model  predictions  directly  to  measured  re¬ 
ceiver  characteristics. 
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•4 


^1.0  FORMULATION  OF  ELECTROMAGNETIC  SCATTERING  BY 
IMPERFECTLY  CONDUCTING  THIN  FLAT  SLABS 


4 . 1  INTRODUCTION 

Predictions  of  ILS  signal  derogation  based  on  the  assumption 
of  perfectly  reflecting  scatterers  have  proven  moderately  success¬ 
ful  in  a  variety  of  typical  problem  situations.  Yet  the  validity 
of  the  perfect  reflectivity  assumption  is  open  to  question  in 
situations  involving  sharp  angles  of  incidence  on  non-metallic 
obstructions. 

To  extend  our  capability  for  predicting  reflection  from  such 
structures,  we  develop  here  an  approximate  model  of  scattering 
applicable  to  flat  walls  or  platforms  of  arbitrary  homogeneous 
electrical  properties.  The  formalism  is  general  with  respect  to 
orientation  of  scattering  surface  and  is  therefore  suitable  for 
treating  properly  the  polarization  of  scattered  radiation.  How¬ 
ever,  to  avoid  difficult  questions  about  field  boundary  conditions 
on  end  faces  and  at  face  edges,  attention  is  restricted  to  thin 
flat  plane  parallel  slabs.  For  complex  structures  which  may  con¬ 
sist  of  a  number  of  thin  slabs  in  combination,  the  problem  of 
mu It isle  reflections  remains  inadecuateiv  treated. 


4.2  DEVELOPMENT  OF  GREEN'S  INTEGRAL 

Ke  begin  with  the  vector  Green's  equation  for  the  scattered 
electric  field  analocous  to  Equation  (2.31): 


-  i  f  r  -  - 

Egir')  =  -  vr  J  ds jissu  (nxK)  v 

s 

-  (nxE)  x  Tv  +  (n'E)T.J 


(4.13 


where  again  E  and  H  are  total  electric  and  magnetic  fields  on  a 
closed  surface  S  enclosing  the  scatterer,  n  is  the  inward  normal 
to  the  surf-  e,  and 


=  exp  {ik|r ’-rsS )/|r’-rsi  . 


(4.2) 


In  the  treatment  of  scattering  by  perfect  conductors,  the  toms 
in  Ecuation  {4.1)  involvinc  n  x  E  and  n-E  are  discarded  on  the 


crounus  tnat 


r'anisnes  at  tne  sort  ace  anc 


:ives  on iv  a 


longitudinal  field.  Bv  aooeal  to  the  Kir  chef  f  aooroxims 


total  field  K  and  H  are  assumed  to  be  zero  on  the  shadowed  portion 
of  the  surface.  Consequently,  for  conducting  surfaces,  the  scat¬ 
tered  field  is  given  by  an  integral  of  the  ’-current  density”  n  x  K 
over  the  i nominated  surface.  In  the  analogous  treatment  of 
scattering  by  a  dielectric  or  poorly  conducting  surface,  the  inte¬ 
grand  term  (n  x  E)  x  ?v  of  Equation  (4.1)  must  be  retained  and  the 
integration  extended  to  the  back  side  of  the  object: . 

Over  the  illuminated  portion  of  the  surface  the  electric 

field  appearing  in  the  integrand  may  be  identified  as  the  sun  of 

the  incident  field  E-  and  a  reflected  field  Er ;  over  the  shadowed 

portion  the  combined  incident  and  scattered  field  may  be  termed 

the  transmitted  field  Efc.  Since,  in  General,  it  is  virtually 

*-  - 

impossible  to  analytically  determine  E  and  over  the  complete 
surface,  it  might  supposed  that  this  approach  would  be  of 
dubious  utility  in  a  practical  scattering  problem.  However,  w-e 
prcceed  on  the  conjecture  that  perhaps  many  structural  elements 
of  ILS  signal  interfering,  buildings,  hangars,  etc-,  i.e.  rein¬ 
forced  partitions,  may  be  reasonably  modeled  as  thin  plane  parallel 
slabs  with  homogeneous  electrical  properties.  This  seems  to  be  the 
logical  extension  of  our  model  for  conduce ing  surfaces. 

Tne  problem  considered  is  that  of  estimating  the  scattered 
field  at  a  point  r '  remote  frees  a  thin  plane  parallel  slab  of 
finite  dimensions  placed  with  arbitrary  orientation  in  an  incident 
field,  as  shown  in  Figure  4.1.  The  illuminated  face  of  the  slab 
is  designated  S, ,  the  opposite  face  S, .  {The  peripheral  surfaces 
of  the  slab  are  neglected  in  the  ensuing  calculations  and  presumed 


to  contribute  neciicxsie  scattering. 


mis  is  tne  most  serious 


deficiency  of  the  proposed  approximation.)  Positions  and  direct 
are  referred  to  an  origin  O,  conveniently  chosen  in  the  vicinity 


ections 


of  the  scatterer.  The  vector  from  this  origin  to  a  point  ?  on  the 
front  face  S,  is  denoted  by  o;  the  vector  to  the  corresponding 
point  on  the  back  face  designated  as  ?*  is  c  -  On,  where  n  is  the 
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has  been  suggested 
calculated  field  a: 
in  a  far  field  aoo: 


discontinuities  at  the  eages  of  S.  ar.c  S?. 


xes  can  se  cosnnersatea  tor 


rocucine  aiec- 


tharga  distributions  on  the  boundary  contours  as 
by  Stratton.*  The  resulting  corrections  to  the 
re  longitudinal  and  accordingly  -ay  be  neglected 
-oxidation .  Thus,  to  obtain  the  proper  far 
:he  scattered  field  due  to  the  cent:  ir-’fiens  of 


field  estimate  of 
surfaces  S.  and  S~  to  the  integral  {4.1;,  ve  ? i: 

:  ieic  given  by - 


:r.e  transverse 
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conditions  for  plane  waves: 

E  =  E.  +  E 

i  i 

g  =  Si  *  K  O"  S, 
K  =  H.  on  S_ 


(4.7) 


In  this  “locally  plane  wave  approximation* ,  the  incident  radiation 
from  a  distant  localized  source  is  aoorcxinated  near  the  coint  ? 


Oil  S 


a  crane  wave  orooacatmc  in  the  direction  or  a  unit  vector 


r. . ;  the  reflected  wave  emanating  free  ?  propagates  in  a  direction 
given  by  and  a  transmitted  wave  emanating  from  face  also 
propagates  in  the  incident  direction  n^.  The  relationships  between 
incident,  reflected,  and  transmitted  fields  will  be  expressed  in 
terms  of  certain  dyadic  forms  which  hold  for  the  boundary  value  pro¬ 
blem  of  the  plane  wave  incident  on  an  infinite  plane  slab.  If 

E  (c)  is  the  incident  electric  field  at  point  ?,  the  local  plane 
o 

wave  aoproximate  fields  are 


i  (kn-  *  Cr-oi  -  sit} 


Si  CD  =  EctDe 


.  _  _  _  i(kn  -(r-o)  -  -t) 

E^Cr)  =  ?-  -  E„(De 


(4.3) 

(4-9) 


•  (r- 


r*&)  -  ut> 


:)  =  T  -  E  (c>e 
o 


(4.10) 


where  F,  and  T  are,  as  yet,  undefined  reflection  and  transmission 
dyadics.  In  the  following ,  the  periodic  time  dependence  vd.ll  be 
suooressed,  but  is  still  understood. 


Taking  advantage  of  Equations  (4.4)  th 
the  relations: 


uch  (4.7)  and  also 


imt  H  =  ikn  x  E 
r  r  z 
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we  are  able  to  reduce  Equation  (4-3)  to  a  single  surface  integral: 


Es(r*) 


-  ik  i  ^  )  ► 

rM  =  1  ^  > 


r  x 


.  f 

r*  x  y  ds 


n  x  (a;  x  E 

r  - 


i  +  Sr  x  Er)  -  ^n  x  (2.+Er)^x  r‘j< 


iA(r*|p) 


£nixEt)  -  (nx£t)  x  r*j  el4{r‘ ^ 0+Dfi) j  I  (4-12) 


where  the  transmitted  field  associated  with  the  position  vector 
o  is  evaluated  at  the  point  o  +  Bn  of  surface  -  This  result  can 
be  expressed  algebraically  as  an  operator  acting  on  the  incident 
electric  field,  if  we  now  introduce  the  dyadic  relations  (4-3) 
thrones  (4-10): 


i s£->  -  ii|irO  x  !;.  xf 


ds 


+  £nr*+a,.n-S-  (r'-n^J  ”1  V  el£*r'  iP> 

-[SJ-«.;-n-(r-«i)l]Teii<;,l»4l,ai)  .  eo(3)  (4.13! 

The  dyadic  expressions  here  have  the  usual  linear  properties;  for 
example: 

(nr'-sn-n)  -  X  =  n(f*  -X)  *  n.(n-X), 


etc.,  -here  x  is  an  arbitrary  vector;  ti 
noted  by  I: 


i-  -  (4.1 4} 

icentity  dyadic  is  de- 


I  -  X  =  X 

T  •  "s 


(4.15) 


To  reduce  Equation  (4-13)  to  its  most  convenient  fora,  we  assume 
that  the  largest  slab  dimension  is  much  smaller  than  the  distances 
to  either  the  radiation  source  or  the  point  of  detection. 
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This  condition  must  hold  if  the  far  field  approx imat ions  being 
used  are  to  be  valid.  This  being  the  case,  we  will  neglect  the 
variation  of  n-,  n  ,  acid  r '  over  the  slab  face  S.  and  set  each  of 

X  T  X 

these  vectors  to  its  respective  direction  at  the  fixed  reference 
point  O.  we  shall  also  expand  the  phase  correction  i  (r 1  j  c+Dn)  and 
k jcp  only  the  love s t  order  terns  in  the  dimension  D,  i.e.: 


fiU’jo+Dn)  =  £{r*  jo)  -  kD(r*  *n) 
fhis  siaolif ication  oeraits  us  to  write  (4.13)  as: 


EsCr*>  = 


vv  )  J(nitr‘tn)  -/^dsHo(a] 


M( r*  !o) 


(4.16) 


(4.17) 


&  (L.r' ,n)  =  r*  x  r*  x  frnr'+nJi  -  a- (r'+nJI^ImR-^T^ 
+  2(n-n.  i(I-nn)  R^J 


with 


and  where  we  have  used 


-ikD(r  *  *n) 

=  e 


n  =  (I-2nn)  -  n. 
r  x 


(4.13) 


(4.19) 


(4.20) 


The  application  of  the  foregoing  results  to  scattering  by  a 
serfectlv  conducting  slab  can  be  made  bv  cutting 


R  =  2nn  -  I 


T  =  0. 


(4.21) 


Equation  (4.13)  then  reduces  to 


V;'>  ■ 


f2?'  ^[(»-ioI;))6.  -  (a-fii,Eo(;i].“‘?,i^ 


(4.22) 


which  result  can  be  obtained  by  straightforward  manipulations 


«irectlv  from  (4.3). 
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To  extend  Equation  (4.13)  to  the  general  case,  we  must  now 
make  a  detailed  analysis  of  the  slab  boundary  value  problea  to 
determine  the  dyadics  R  and  T. 


4.3  REFLECTION  AND  TRANSMISSION  0?  PLANE  HAVES 
BY  AN  INFINITE  PLANE  PARALLEL  SLAB 


In  order  to  apply  the  local  plane  wave  approximation  in  our 
estimates  of  scattering  by  thin  slabs,  we  solve  the  boundary  value 
problem  for  a  plane  wave  incident  on  an  infinite  plane  parallel 
slab  of  thickness  D,  with  homogeneous  electrical  properties  e,  u 
and  a  which  separates  two  regions  of  vacuum.  Because  of  the  plane 
symmetry,  the  problea  is  two-dimensional ,  with  the  incident,  re¬ 
flected,  and  transmitted  wave  normals  at  both  slab  surfaces  all 
parallel  to  a  common  plane  of  incidence.  The  incident  wave  of 
arbitrary  polarization  can  most  easily  be  treated  by  decomposition 
into  an  "E  "  ccrusonent,  with  the  electric  vector  normal  to  the 

4.  * 

plane  of  incidence,  and  an  "E i I "  component ,  with  the  magnetic  vec¬ 
tor  normal  to  the  plane  of  incidence.  '  We  treat  the  ea  case  first. 

The  several  secondary  waves  generated  by  an  electromagnetic 
wave  incident  on  a  slab  boundary  are  illustrated  in  Figure  4.2. 
On  the  incident  side  of  boundary  ,  w*e  have  the  reflected  wave  E 
in  addition  to  the  incident  wave  £..  In  the  slab  medium,  we  have 

_  i  ^ 

a  "forward  wave"  E,  and  a  "backward  wave"  E,  .  Outside  the  slab  on 

f  b 

the  transmission  side  of  S2»  we  have  E^ .  All  of  these  fields  are 

normal  to  the  plane  of  incidence,  and  we  nay  drop  vector  notation 

for  these  components.  We  denote  the  value  of  a  field  at  a  slab 

boundary  by  a  subscript  1  or  2;  the  wave  normals  are  denoted  by  the 

vectors  n.,  n  ,  n_,  n.  .  and  ri  (=  n, ) ,  as  shown  in  Ficure  4.2.  A 
i  r  z  o'  t  i 

plane  wave  in  the  slab  has  a  complex  propagation  vector  £.  The 

A  A 

real  part  has  the  direction  of  either  r,f  or  n,^  and  determines  the 
surfaces  of  constant  phase.  The  imaginary  part  is  parallel  to  the 
boundary  normal  n  and  determines  the  surfaces  of  constant  amplitude. 
Outside  the  slab,  the  propagation  vectors  krL  and  krv  are  real, 
with  k  =  23/1- 

Following  Stratton, x  the  boundary  conditions  at  surfaces  S, 

and  £,  for  the  E,  polarization  can  be  written: 

2  * 
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J 


Ei  +  Er  -  Ef]_  +  Eb2 


(4.23) 


i  (ln-i:i)Ei  ♦  (K-Er)Er)  4[(S^f1Efl  +  (S-?b)Ebl]  i4‘24> 


Ef2  +  Eb2  "  Et 


(4.25) 


i  [<n.K£)Ef2  *  (fi-ICb)Eb2]  =  ^iEt 


(4.26) 


Two  additional  conditions  on  the  six  undetermined  variables  are 
needed;  they  may  be  taken  to  be: 


where  we  have  used 


i  (n-Kf  )D 

Ef 2  =  e  Efl 


=  6E, 


Ebl  ”  B  Eb2 


ri*Kf  =  -  n.Kb 


(4.27) 


(4.28) 


(4.29) 


which  generally  holds  for  plane  wave  reflection.  By  making  use  of 
the  similar  relation 


n-k.  =  k  n*n. 

l  i 


=  -  k  n*n 


(4.30) 


we  can  reduce  Equations  (4.23)  through  (4.28)  to  the  solution  for 
the  normal  electric  sector  components: 


t  =  - 


(1-G2) (1~V^) 


[{1,v2  -  b2u-Ye>2]  1 

-  Vi 


(4.31) 


[(.vte)2  -  B2(1-ye)2] 


21  Ei 


(4.32) 


=  TE  Ei' 


iittSiiilii.Ni  wuiU'iw  m  Ml  Mill  to  lljfldiftlidrtlAS1  WS  #1  <f»gWK  fgflwj  >  npfa,*  jjyjjjk  ■ « AMteifawstf w 


(4.33) 


We 


■y  =  y  (fi;fy _ 

E  v0(n-tf) 

->  A 

wish  to  calculate  the  quantity  Kf*n.  We  have  the  expansion 


Kf  =  n{n-^f)  -  n  x  (nx"£f )  . 


(4.34) 


We  also  have  the  boundary  condition 

n  x  =  n  x 

where  k.  is  the  incident  wave  vector. 

i 


Then 


=  (n*rf)2  +  j^nx  (nxki)  j  *  |fi  x  (nxk^J 


«  (n*Kf)2  +  (nxki)  *(nxk±) 
=  (n-Kf)2  +  k?  sin2  9± 


(4.35) 


(4.36) 


But  we  have  also,  as  shown  by  Stratton, 

2  2 

=  as  ye  +  l  onio. 
From  the  above  equations,  it  follows  that 
(n  •  J^)  =  f 


2  ,2  .  2  . 

as  wt  -  k  sm  9^  +  l  as ua 


(4.37) 


(4.38) 


hence 


ye  =  'i 


k  cos  6 . 

i 


6  =  e 


°  ^  a? v  (e+ic/cu)  -  k^  sin2 

r  35  $  O 

iD  (e+in/sa)  -  sir.'4'  5^ 


(4.39) 


(4.40) 


For  the  Ei ,  polarization  case,  the  magnetic  vector  is  normal 
to  the  plane  of  incidence.  Evaluation  of  the  transmitted  and  re¬ 
flected  waves  follows  from  the  set  of  boundary  conditions  similar 
to  (4.23)  through  (4.28): 

Hi  +  Hr  =  Hfi  +  Hbi  s-41> 
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(4.42) 


(H.-H 

l 


n-Kf 
e+i a/® 


) 


Hf 2  +  Hb2  =  Ht 

(4.43) 

Kf 

_.-i  /u  {Hf2  Hb2}  c  Ht 

(4.44) 

Hf 2  *  3  B£1 

(4.45) 

Hbl  *  8  “b2 

(4.46) 

We  obtain  results  analogous  to  Equations  (4.31)  and 
with  y  _  replaced  by  y„; 

t*  n 

(4.32)  above 

(e+io/u)  n-k^ 

V'f 

(e+i o/u) 

u  ec  'e 

(e+ic/u)k  cos  6^ 

,  2  2 
£q  Ww*'u  (e+ic/w)  —  k  sin  6^ 

(4.47) 

Hr  =  »H  Hi 

(4.48) 

Ht  "  TH  Hi 

(4.^9) 

(1-yJ) (1-B2) 

rh  “  — — 5 - 5 - 5 - 

(1-B  )  (1+yJ)  +  2(1+B  )Yh 

(4  50) 

4y.,S 

iji  ^  *• » 

H  (1-B2)(1+y2)  +  2(1+62)vh 

(4.51) 

To  use  (4.31),  (4.32),  (4.48),  and  (4.49)  in  the  general  case, 
we  must  decompose  the  incident,  reflected  and  transmitted  electric 
vectors  into  polarization  components  normal  and  parallel  to  the 
plane  of  incidence.  For  example,  the  component  of  an  incident  wave 
perpendicular  to  i^oth  the  boundary  normal  li  and  the  wave  normal 
may  be  expressed  by  the  equation: 


Ei^  =  — 2 —  (nxrK)  (nxiK-E.^ 


(4.52) 


sin  s ^ 


In  anticipation  of  dyadic  relations  to  be  developed  shortly,  it  is 
appropriate  to  introduce  a  special  dual  set  of  basis  vectors  to 
express  such  polarization  relations.  Consider  the  linearly  inde¬ 
pendent  basis: 


el  sinS . 


n. 

-  x 

e2  ~  sinB^ 


n  x  n- 

e  =  — . — . — — 
3  sins^ 


(4.53) 


(see  Figures  4-3A  and  4.3B). 

These  vectors  are  dual  to  the  set 


-1 

e  = 


uTsT  (n“(n'ni)Ri) 

.  — fn . - (n-n . )n\ 

V.  x  i  / 


in  the  sense  that: 


1  -  ~  ~  \ 
UnsT  (nxV 


**■  -*-1 

e  -  *e  =  o .  . 
i  13 


(4.54) 


(4.55) 


The  set  {e-5}  is  the  appropriate  basis  for  representing  the  incident 


electric  field  on  a  medium  interface  normal  to  n.  If  we  make 


E±  =  (Ei)j  eJ 


(4.56) 


where  the  repeated  Roman  index  here  and  in  the  following  will  imply 


a  summation 


3 

I 

3=1  , 


then: 


Figure  4.3A.  Geometry  for  the  Vectors  ej_f  &2  and  e3*  x'»y**z 

Form  a  Right-Handed  Triad,  je.'  and  y*  Lie  on  the 
Surface  of  the  Slab,  z*  ,  n^  and  fi  Lie  on  the 
Plane  of  Incidence.  Then  £3  is  in  the  Direction 
of  x*  and  is  £  to  the  Plane  of  Incidence.  (z*  i 
Not  Necessarily  in  the' Vertical  Direction.) 


/ 


+ 1  +2  +3  ■  a 

Figure  4.3B.  Geometry  for  the  Vectors  e  ,  e  and  e_^.  z  ,  r.j_ 
and  n  Lie  on  the  P^ane  of  Incidence,  e3  is  1  to 
this  Plane,  While  e1  is  Parallel  to  it  and  is  1 


(4.57) 


(E.) .  =  e.  *  E. 

ID  D  i 

Upor.  comparison  with  Equation  (4.52),  we  see  that  eJ  (e^  *E^)  is  the 
perpendicular  component  E.  ,  that  (e2'S.)  is  always  zero  for  a 
transverse  wave,  and  the  e  (e^*E^)  is  the  parallel  polarization 
component  of  the  incident  wave. 

The  reflected  and  transmitted  electric  vectors  for  each  polar¬ 
ization  are  now  easily  calculated.  First,  in  the  perpendicular 
oolarization  case  we  have: 


ErJ_  =  *£  e3(i3-Ei) 


E  =  T  e3  -=■ 


'Ov8i) 


(4.58) 

(4.59) 


In  the  parallel  polarization  case,  we  have 

k 
'o 


E  I  i  - - —  n  x  H 

r‘  ■  ioe_  r  rj 


Ru  n  x  H. 
H  ut  r  x. 
o 


-1  ~  - 


=  -  Rjj  nr  x  (n^xe  )  (e^E^) 
By  a  similar  analysis,  we  find: 


Efc  |  |  =  “  th  "i  x  (n^xe1) (e, -E . ) 


i  1 


(4.60) 


(4.61) 


Since  the  sum  of  the  two  polarization  cases  gives  the  total  re¬ 
flected  and  transmitted  fields,  we  are  now  able  to  write  down  the 
+ — ► 

dvadics  R  and  T: 


R  =  -  Rh  nr  x  (nixe1)e1  +  Rg  e3e3 


T  =  -  T  n.  x  (n.xe1 )e.  +  T  e3e, 

n  1  X  X  £>  j 


(4.62) 


(4.63) 


Here,  nr  is  the  wave  normal  of  the  reflected  wave  generated  by  the 
given  incident  wave  in  the  local  pla..e  wave  approximation,  i.e.: 

n  =  n.  -  2(n-n.)n  (4.64) 

r  i  i 
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It  can  be  checked  again  that,  in  the  case  of  a  perfect  conductor 


R_  —  1 
R..  -  ♦  1 

n 

[-S2?!  *  (4.72) 

as  they  should. 

4.4  SCATTERING  FROM  A  SLANTED  RECTANGULAR  CONDUCTIVE  SLAB 

In  this  section  we  will  calculate  the  far  field  scattering 
from  a  slanted  rectangular  conductive  slab,  the  to?  and  hot ton 
edges  being  parallel  to  an  infinitely  conductive  flat  ground  plane. 
The  general  method  given  in  Sections  4.1  and  4.2  will  be  applied 
to  this  problem.  Figures  4.4A,B  and  C  show  the  geometry  for  this 
problem,  n^,  the  incident  wave  direction,  is  approximated  by  R', 
the  unit  vector  from  the  antenna  to  the  reference  point  0  on  the 
midpoint  of  the  slab's  lower  edge. 

rL  =  R*  =  i  cosS^  cosc  +  j  cosc^  sins  +  £  sinSj  (4.73) 

where 

sin&1  =  (M-Ha)/R' 
cos “i  =  ddi/r’ 


R'  = 


[D?l 


(m-ha)2]1/2 


The  incident  electric  field  vectors  (with  the  inclusion  of 
ground  reflection)  for  small  elevation  angles  is 


E  (o)  =  U  E  f(u) 
o  o 


sikR 

R 


1-e 


ik(2z'VDP) 


(4.74) 


(Compare  Sc.  (2.58).)  Here,  the  angle 


v  =  tan-1  (Vj^/Xj^), 


x,  and  y,  being  the  horizontal  coordinates  of  the  reference  point 
0  on  the  slab  with  the  horizontal  edge  of  length  L.  Ke  make  the 
approximation  and  again  from  the  expansion 
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Figure  4.4A-  General  View 

xrnimd  plane . 


R  =  R*  + 


R**  0  +  [o2  -  (R*-p)2] 


/2R*  -f- . 


(4.75) 


we  obtain  (supposing  than  ko2/2R*<<l)  the  Fraunhof ffer  approximation 
for  the  incident  field 


ikR  ikR* 
e  -  e 


"5t- 


ik (R* -s) 


(4.76) 


From  Equation  (4.22)  we  have  the  scattered  field  from  a  conductive 
flat  surface 


;  _  u  e 

's  2- 


ikr 1 


r*  x  |^(nin-  )  -  (n-n^  J dsEQ  (o) 

S, 


x  e 


id(r*/o) 


\  « 


(4.77) 


where  is  the  illuminated  front  surface,  and  where 
=  ~  kr»-o  +  k[o2  -  (R-S)2]/2r*  + . 


(4.78) 


In  the  calculation  in  this  section  we  assume  the  Fraunhofer 
field  condition  to  hold  for  the  scattered  field. 


.  2  _  2 

^r=  «  1 

2r '  <r 


(4.79) 


so  tnat 


Then, 


where 


&(r’/o)  =  ~  k  r‘ 


Fs  =  r'  xpitn-O)  -  U<n.5.)j}lb 


ikZ  f(£> 


I  =  — —  exo  ^ik(r'+R’)  I  /r*R* 
a  t»  J  c 


(4.80) 


(4.81) 


(4.82) 


ana 


= J  dn  dC  ^exp  jik  (R*  *o)j  ^1-exp  (ik2z *  exp  {-ikr  *  *o)| 


(4.83) 
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with 


sine2  =  (z2-M)/r’  (4.915 

Also  la  Equation  (4. SI)  we  nave 

n^n-U)  -  (n-n^U  =  n  x  (n^xU)  (4.92) 

The  induced  voltage  at  the  aircraft  antenna  is 

V  =  cQ (kxr * )  -  Es  (4.92) 

V  =  cq  | Coer' )  *  r*  x  J^r*  x  (n  x  (n^O))]}  Ib  (4.94) 


where  I.  is  civen  bv  Eauations  (4.82)  and  (4.87).  Ke  have. 

O' 

finally, 

V/c0Ifa  =  |[{kxr*)  x  r*J  x  r*}  - 1  n  x  (n±xU)] 

=  -  r*  x  k  ♦  ^(n^xU)  x  nj 
=  -  r*  -  cosS^  +  sinS^  (kxU)k-n 

-  ccsSj^  k(k*n)|  (4.95) 

Assuming  the  unit  vector  U  to  be  horizontally  polarized 

U  =  -  (kxn^)/cos5^  (4.96) 

we  obtain  fer  the  voltage  function 


V/c  I-  =  cose.,  [cosr,  sine  cos 5.  +  sins,  cose  cos  ( 
on  1  x  1 

Si+Yi)]  (4.97) 

In  order  to  obtain  the  voltage  due  to  the  direct 
the  wall,  we  set  (see  Figure  4.5): 

scattering 

sinc2  =  (z2~H)/r ’ 

(4.93) 

«•  -  [■&  *  <’2->2]1/2 

(4.99) 

CCS  c2  =  t»a2/r* 

(4.100) 

SS 

REFERENCE  PT.  ON  SLAB 


OBSERVATION 
POINT  AT 
AIRCRAFT 


Figure  4,5.  Geometry  for  Calculating  the  Scattered  Field  from 
the  Ground  Image  of  the  Slab 
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In  order  to  obtain  the  scattered  electric  field  from  the 
image  of  the  wall,  we  first  formally  compute  the  scattered  field 
from  the  wall  to  a  point  below  the  ground  plane,  which  is  the  same 
distance  as  the  observation  pcilit  above  the  ground  plane,  and  then 
"reverse"  the  horizontal  component.  Since  in  the  voltage  equation. 


v  =  cQ(kxr')  *  Es, 


(4.101) 


k  x  r*  is  a  horizontal  vector,  after  formally  calculating 
the  field  at  the  point  below  the  ground  plane,  one  can  simply 
multiply  this  quantity  by  -1,  since  the  vertical  component  is  can¬ 
celled.  Figure  4.5  shows  the  geometry  for  this  image  calculation. 
The  image  field  at  r"  is  obtained  by  setting 


sine2  =  -  (M+|z2! )/r' • 


where 


r"  =  d£2  +  (M+ !  z2  I ) 2  1/2 


cose2  =  Dp2/r" 


(4.102) 


(4.103) 


(4.104) 


4.5  SCATTERING  FROM  A  SLANTED  RECTANGULAR  DIELECTRIC  SLAB 

The  geometry  for  the  scattering  is  show  in  Figures  4.4A,  B 
and  C.  Many  of  the  quantities  calculated  in  the  previous  section 
(4.4)  for  the  conductive  case  apply  to  the  -eneral  case  treated 
here. 

The  scattered  electric  field  is  given  by  Equation  (4.17), 
which  can  be  written  as 

•*— * 

£ ,  , ,  ik  exD  ( ikr ' )  cz ,  ~  ~ ~ . 

E  (r  )  =  — *— p - J  (n±  ,r  ,n)  - 


*  f  ds  EQ(p)  exp  |i&  (r  '/p)J 


(4.105) 


■ W  V  J****~***w&%tt  . . . . » *********  imm » . . . . . . >)l|lll||^^ 


With  the  same  far  field  approximations  as  in  the  previous  section, 
we  have  for  the  incident  field  on  the  slab, 

£  (o)  =  F  U  (4.106) 

o  n 


where 


U  =  (e. -U)e '  +  (e,-U)e3 

i.  J 


and 


eikR'  ik£*.' 

Fn  S  Eo  f{*>  V“  e 


2ikz 'H./D 
.  A  p 

1-e  r 


Then,  defining 


where 


with 


we  have 


ra  '  /ds  5o(0)  eU“”  *  «  xi 


iA(p) 


rb  ‘  /  ds  Fn  e 


h(ryp)  =  -  k  r*  *  p 


ikR* 

*b  =  Eo  f  W  V-  Jc 


(4.10-’) 

(4.1C8) 

(4.109) 

(4.110) 

(4.111) 

(4.112) 


where  Ic  is  given  by  Equation  (4.87)  in  Section  4.4.  Combining  the 
previous  equations,  we  obtain  for  the  sca'ttered  field 

Es  =  F2  (5-U)  (4.113) 

where 

-  ik(r'+R’) 

F2  -  k  Eo  f  <*>  S-R1 -  1c  (4*114) 

The  induced  voltage  at  the  aircraft  antenna  is 

V  =  C  (kxr ' ) -2 
o  s 

*4— ► 

-  Co  F2(kxr') . (7-U)  (4.115) 
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Referring  to  Equation  (4-71)  we  haves 
y  -u  =  r*  x|r*  x  ^(iixe3)  (-l+OjjTn-Rj.) 

+  (n*niX-i+aDTK+RH)  (r’x^3))  (1,-U) 

+  ((r -xe2)  (-l+Vrt) 

+  {n*ni)e3  (-1+a^+i^))  (e3-U)]| 


(4.116) 


If  we  denote  the  quantity  enclosed  by  the  square  brackets  t  j 
by  A,  the  voltage  function  may  be  written  as 


V/CqF2  =  (kxr ' )  •  r*  x  (r'xA) 

This  can  be  written  as  a  sum  of  four  parts, 

&7 =  (7^)1 +  fe), +  +  (277), 


We  have 


e3  sinS^ 


e^  =  n/sin8^ 


n  =  £  sine  sina  +  3  sine  cosa  -  k  cose 


(4.117) 


(4.118) 


(4.119) 


(4.120) 


(4.121) 


1*2  -  1  cose2  sin(B^  -  a)  -  j  cose2  cos  (3^  -a)  (4.122) 


+  k  sine- 


After  some  algebra  we  obtain 


/  XT  \  A  A>  A  A  A  AAA 

fc~F~j  =  “  (n*ni}  k  *(r'xn>  +  k-(r'xni) 


-  ~  ,  1  {-1+aDTH_RH)(el*U) 

*  (r'xn)  +  k-  (r  'xn. }  - -  ■■■-5  ■ - - - 


sin9 . 
1 


(4.123) 


lilllillll'IlllNilliliKItliillii  Niili.llilliU  1 1 
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where 


k'(r‘xn)  =  cosc2  sins  sinS., 


(ii*R.)  =  sine  cos 5 .  cosy-]  -  cose  sinS, 


k-  (r'xn^)  =  cose2  cosS-j^  sin^  +  y1) 


(4-124) 


(4.125) 


(4.126) 


sin'3  ^  =  j  nxf.±  J 


(.2.2  2  2  2  2  2 
=  <sm  3^  sxn  e  +  cos  cos  e  +  cos  9^  sin  e  sin  y^ 


+  2  sinS^  cos9^  sine  cose  cosy^ 


(4.127) 


Also  we  obtain 


(V  \  #Ai  A  i  ^  A.  *  A  A  A  1 

CT^)  ~  ~  tr  niJ  k*(r  ~  (r’*n)  k.  (r’xn^ 
X  [("*"i}  t~1+aDTH+RH)  (®l'S)]  / Sin8i 


where 


r**n^  =  -cose,  cos  S^cos (Bj+y^)  +  sine2  sin6^ 
r'*n  =  -cose2  s^ne  cosB^  -  sine2  cose 

and 

-  (£.£.)  | (r’-n) (r'*k)  -  (n-k)}] 
x  [{-1+aDTE'RE)  /sin6i 


(4.128) 


(4.129) 


(4.130) 


(4.131) 
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-*  i  r- 

e,  =  — = — r—  ji(sin6,  sine  coso+  cosS,  cose  sm?) 
3  smc .  l  1  1 


-  j (sinS^  sine  sina  +  cos9^  cose  cos^ 


-  k(cosS1  sine  siny^)J 


(4.139) 


Accordingly,  when  the  components  of  u  are  determined  the  products 
<^-u  and  e3*u  can  be  computed  for  the  general  case  of  a  non-hori- 
zontally  polarized  electric  field  incident  on  the  slanted  di¬ 
electric  slab. 


APPENDIX  A-  THE  FRESNEL  APPROXIMATION 


I -  INTRODUCTION 

Equations  (2.47)  and  (2.48)  of  this  report  describe  the  scat¬ 
tered  magnetic  and  electric  fields  in  the  Fraunhofer  zone  of  a 
perfectly  conducting  scatterer.  These  equations  represent  the 
leading  terms  in  the  asymptotic  expansions  of  Equations  (2.39)  and 
(2.42)  for  large  values  of  the  distance  between  the  scatterer  and 
the  receiver.  Let  D  denote  some  characteristic  linear  dimension 
of  the  scatterer,  and  let  r*  denote  the  distance  from  some  point 
on  the  surface  of  the  scatterer  to  the  receiver.  The  assumptions 
involved  in  the  Fraunhofer  aonroxination  can  be  summarized  as 


follows. 


r*  >>  E» 


D  /2r ’  «  X 


(A.l) 


(A. 2) 


where  X  is  the  wavelength  of  the  incident  radiation.  In  this 
appendix,  condition  (A. 2)  will  be  relaxed,  and  Fresnel  zone  ap¬ 
proximations  for  the  scattered  fields  will  be  developed  and  then 
applied  to  the  problem  of  localizer  signal  scattering  by  a  flat, 
vertical  wall.  The  motivation  for  this  new  approximation  is  the 
ongoing  construction  of  very  tail  buildings  (motels,  747  hangars, 
etc.)  near  airport  runways.  For  such  structures,  condition  (A. 2) 
is  often  violated,  at  least  over  portions  of  the  aircraft  approach 
paths. 

II.  ANALYSIS 

The  asymptotic  expansions  of  Equations  (2. 39)  and  (2.42)  for 
large  values  of  r*  were  performed  upon  the  phases  and  amplitudes 
of  the  two-point  Green’s  function  y(r*,r)  and  its  derivatives 
which  appear  in  those  equations.  The  function  ?  is  given  by 


f(r’.r)  = 


iklr’-r! 


jr'-rl 


(A. 3) 


A-i 


where  k=2 -/•  is  the  wave  number  of  the  Incident  radiation,  and  r* 
and  r  denote  respectively,  the  position  vector  of  the  receiver  and 
the  position  vector  of  a  point  on  the  illuminated  surface  of  the 
scatterer  relative  to  an  origin  of  coordinates  0  which,  for  con¬ 
venience,  ve  locate  on  the  surface  of  the  scatterer.  Let  D  denote 
the  largest  value  assumed  by  r= ! r  j -  D  represents  sene  character¬ 
istic  dimension  of  the  scatterer-  If  r ' >>D  (condition  (A-l)*, 

- .  -1 

tne  anplituae  JLr’-rj  of  y  can  be  approximated  very  accurately 
simply  by  (r*)_i.  More  care  must  be  exercised  in  approximating 
the  phase  of  y  because  of  the  oscillatory  behavior  of  the  complex 
exponential-  Expanding  i r  * -r  j  to  terms  of  second  order  in  r,  we 


(A.4) 


where  r*=r’/r*  is  a  unit  vector  in  the  direction  of  r'.  If 
■> 

D  /2r*<<\  (condition  f A - 2 • ) ,  we  can  clearly  neglect  the  quadratic 
terns  in  (A.4)  and  use  the  first  two  terns  in  approximating  tne 
phase  of  y  for  large  values  of  r*.  This  “linear"  representation 
for  the  phase  of  ?  is  just  the  Fraunhofer  approximation  used  in 
Section  2-0  of  this  report,  if  D2/2r*  is  not  snail  compared  to  *, 
the  quadratic  terms  in  (A.4)  cannot  be  ignored,  and  ve  arrive  at 
the  following  asymptotic  representation  for  ?r 


Equation  (A-5)  is  the  asymtotic  representation  of  the  Green's 
function  =  in  the  Fresnel  zone  of  the  scatterer-  It  should  be 
emphasized  that  in  deriving  Equation  CA.5),  we  have  assumed  that 
the  phase  or  the  Green's  function  can  be  accurately  represented  by- 
expanding  |r*-rjt  to  terms  of  only  second  order  in  r.  "that  is  to 
say,  it  has  been  assumed  that  all  higher  order  terns  are  snail 
compared  to  a  wave length. 

when  the  asymptotic  expression  for?  given  in  Equation  (A. 51 
together  with  the  corresponding  asymptotic  expressions  for  the 
various  derivatives  off  are  substituted  into  Equations  (2-37}  and 
(2-39),  the  fol lowing  asymptotic  expressions  for  the  scattered 


A— 2 


fields  ia  the  Fresnel  zone  of  a  cor  tract  me:  scatterer  are  obtained 


HsCr*) 


=  -i | sill j-,/  <£a.)«'ik,,?’*J,as . 

2T  1 


CA.6] 


M'>  =£(?) 


„xl/2  eikr* 


r  x 


rm  xj  (nxH.)e  ~c:  'r'ds 

»—  -A. 


SM  -  / 


! (r* ,r)  =  r*  -r  - 


ifsere 

!  -  ..  .. 
P-  "  2r"‘  I  *-‘-8 

j 

la  Equations  (A. 6)  and  (A. 7),  S+  denotes  the  side  cf  tlx:  surface  : 
the  scatterer  directly  exposed  to  the  incident  radiation,  5.  is 
the  incideat  atretic  field  ca  SA  arc  a  is  the  unit  inward  noraal 
to  S+  at  r. 


III.  LOCALIZER  SIG!iAL  SCATCERIK'I  BY  A  FLAT,  VERTICAL  *ALL 
H.  THE  PRESSES,  AFPROXIXATICS 


Is  Section  2.5  the  prohlea  of  localizer  signal  scattering  ir 
a  flat,  vertical  wall  was  analyzed  using  the  Fraunhofer  appre-xi- 
saticc.  This  sane  problem  will  new  be  treated  using  the  Fresnel 
equations  (A. SI  and  CA.7) .  The  geometry  of  the  problets  is  illus¬ 
trated  is  Figure  2.5. 


It  was  shows  is  Suction  2.5  that  the  total  incident  sagnstic 
field  H;  (direct  plus  ground  reflected)  at  as  observation  point 
t*,y,z3  in  the  far  field  of  a  horizontal ly  polarized  localizer 
antenna  locates  at  (0,0, E)  could  he  represented  as  fellows: 


a  — 


(t) 


1/2 


ilR  T  >ikz=/Dj 


E  f  (#) 
o 


(A.®: 


where 


(y/x) 


(A. 12) 


In  Equation  (A. 9) ,  f {: )  is  the  horizontal  antenna  pattern  of  the 
localizer,  Eq  is  an  amplitude,  and  e2  is  a  unit  vector  along  the 
z-axis  which  is  perpendicular  to  the  perfectly  conducting,  flat 
ground  plane  {the  x-y  plane).  Equation  (A. 9)  accurately  describes 
H.  at  observation  points  with  snail  elevation  angles  (z/D  <<1) . 


As  in  Section  2.5,  we  will  use  the  n id point  of  the  base  of  the 
wall  as  an  origin  of  coordinates  for  the  surface  integrals  in  (A. 6) 
and  {A-?}-  The  coordinates  of  this  point  will  be  denoted  by 
ix. j,y,,0).  The  azisuth  angle  of  (x^,y, ,0)  will  be  denoted  by  v 
C.  =  tan  ~ (y./x. ) ) .  »fe  again  assure,  as  in  Section  2.0,  that 

n  i 

O  -  =  %'X.  V- 

pi  T  1  ■  i 

is  such  greater  than  the  length,  L.  of  the  wall,  so  that  f {: )  can 
is  reel  a  red  by  f  ( )  and  D  bv  D  ,  in  (A.  9)  for  all  soints  on  the 

*  _  p  pi 

wall.  The  vector  A  free  the  localizer  to  any  point  P  or.  the  wall 
can  be  represented  as  follows; 


3=3,  +  r 


(A. 13) 


where  a.  =  x,e  v. e  -  He 
x  i.  X  x  v 


r  is  a  vector  in  the  olane  of  the 


wall  drawn  frees  the  point  (x.  ,y_ , 0)  to  the  point  ?.  Assuring  that 

=  =  *  v?  +  E4 

i  ri  -i 

is  rich  larger  than  the  largest  dirxension  of  the  wail,  the  distance 
?.  frees  the  localizer  to  P  can  be  represented  to  terns  of  second 
order  in  r,  as  follows; 


iV.a  ■ 


t  2  -4-  ^ 

"1  *  i 


(A -14) 


where  =  R./E,.  In  Section  2-5,  the  cuadratic  oath  ler.cth  dif- 
ference  terns  in  (A. 14)  were  neglected  {ef.  Eg .  (2.54)) .  For  the 
present  analysis,  these  terns  will  be  retained  for  the  sake  of 
consistency,  since  the  tems  ef  second  order  in  r  in  the  series 
expansion  of  tlx:  Green’s  function  7  are  being  retained  (ef.  Eg  - 
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The  vector  r  from  (Xj^/Yj^/0)  to  the  point  P  on  the  wall  can  be 
represented  as  follows: 

r  =  nn  +  2  e2  (A. 15) 

where  n  is  a  unit  vector  in  the  plane  of  the  wall  and  parallel  to 
the  ground  (Fig.  2.5)  and  n  is  variable  ranging  frees  -L/2  to  +L/2. 
The  variable  2  is  just  the  elevation  of  the  point  P  above  the 
ground  and  varies  from  0  to  h  where  h  is  the  height  of  the  wall. 

In  Section  2.5,  the  following  approximate  expression  for  the  dot 
product  (R^-r)  was  derived: 

A  -+■  zfi 

R.*r  =  n  cos  (S4$)  -  - —  ,  (A. 16) 

pi 

where  3  is  the  angle  between  n  and  the  x-axis  (Fig.  2.5).  Sub¬ 
stituting  (A.  16)  into  (A.  14)  and  approximating  by  in 

the  last  two  terms  of  (A. 14),  we  obtain  the  following  approximate 
expression  for  the  distance  R  from  the  localizer  to  any  point  P 
on  the  wall: 


R  :  Rj  +  11  cos  (e+f)  -  —  sin2{6+$)  +  —  ,  (A. 17) 

pl  pl  pi 

—2 

where  terms  of  order  (Dpl)  and  higher  have  been  dropped.  Sub¬ 
stituting  Equation  (A.  17)  for  R  into  the  complex  exponential  in 
Equation  (A. 9)  and  replacing  th*»  factor  R_1  by  Rj”  ,  we  finally 
obtain  the  following  approximate  expression  for  fL  at  a  point  P  on 
the  surface  of  the  wall: 


=  2i  «*  (i)  Eo  £<« 


ikR, 


ijcn  cos(S+^>  kzE 
e  sin  — — 

DoI 


ik|n2  sin2  (3+t)  +  z2]/2D?1 


(A. IS) 


Let  ix2,y2,z2)  denote  the  coordinates  of  the  receiver.  The 
vector  r '  from  the  base  of  the  wall  to  the  receiver  is  given  by 


r*  =  (x„-x1)ex  +  (v,-y, )  e„  4  z^e 


2  51 '  v  *  2  s 


(A. 19) 


Let  R,  denote  r ’  =  jr*  | 
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(A. 20) 


R2  =f(x2“xi)  +  (y2-yi>  +  2  2 

The  horizontal  distance  from  (x^,y^,0)  to  the  receiver  will  be  de¬ 
noted  hv  D  _  : 

p2 


DP2 =  V(vxi)  +  (y2“yir 


(A. 21) 


Assuming  small  elevation  angles  (Dp2>>z2^ '  was  s^own  Section 
2.5  that  the  dot  product  r'*r  could  be  represented  approximately 


as  follows: 


r ' *r  z  n  cos  (y-6)  + 


(A. 22) 


where  y  is  the  angle  between  the  projection  of  r '  onto  the  ground 
plane  and  the  x-axis  (cf.Fig.  2.5).  Substituting  (A. 22)  into 
Equation  (A. 8),  which  defines  the  phase  function  $(r',r)  appearing 
in  (A. 6)  and  (A. 7),  we  obtain  the  following  approximate  expression 
for  4: 


z  z 

0  (r '  ,r)  =  n  cos(v-6)  +  .= — — 

°p2 


b -  sin2(y-e)  -  Id— 
p2  p  z 


(A. 23) 


where  (r')  ^  ^  has  been  approximated  by  second 

and  third  terms  of  (A. 8)  and  terms  of  order  D^"2  and  higher  have 
been  neglected. 

Equations  (A. 18)  and  (A. 23)  for  and  $(?’ »r)  can  now  be 
substituted  into  Equation  (A. 6)  and  (A. 7)  and  the  integrals 
evaluated.  As  in  Section  2.5,  the  range  of  the  z  integration  is 
taken  to  be  -  h  <  z  <  +  h  in  order  to  account  for  the  image  of  the 
wall  in  the  ground  plane.  The  results  of  these  integrations  for 
the  magnitude  of  the  scattered  electric  field  at  the  receiver  is 
given  below: 


ik (R1+R2  ) 

E  «  E  f  ftp)  §-5 - 

s  ir  o  R1R2 


sin(v-0) 


(A. 24) 


. . .  .l  iliji'ln.  Mill  to 


where 


P  r  ■ 

t  f  -  kHz 

11  ~-l  sinLv 


sin 


kz2z 


L°P2 


ik22  (25-7  *  56-r)  ■  (A-25) 

'  Pl  P2'd2  , 


and 


+L/2 


=  f  eikT1[' 
-L/2 


cos(6+t)  -  cos  (y 


-6)] 


ikn 


sin2  (6+£)_  +  sin2  (y-S  ) 


2D 


Pi 


2D 


p2 


Idn 


(A. 26) 


Unfortunately,  the  integrals  1^  and  I2  cannot  be  evaluated  in  closed 
form.  However,  they  can  be  evaluated  numerically  using  series 
representations  for  the  complex  Fresnel  function  F(x)v 


x  ift2 

F(x)  s  /  e  dt 


'/ 


(A. 27) 


IV.  COMPARISONS  OF  THE  FRAUNHOFER  AND  FRESNEL  FORMULATIONS 

In  the  Fraunhofer  approximation  used  in  Section  2.5,  the 
integrals  1^  and  I2  defined  by  Equations  (A. 25)  and  (A. 26)  assume 
the  following  simple  forms: 

h 


kz2z 

DP2 


q  sin  [j^j  sin 

+W2  r  *| 

j.  _  j  eikn  cos(6+iJ/)  -  cos{y-6)J 

-T./0 


dn 


(A. 28) 


(A. 29) 


Of  course  the  integrals  appearing  in  Equations  (A. 28)  and  (A. 29) 

can  be  explicitly  evaluated  and,  in  fact,  are  proportional  to  the 

sine  functions  appearing  in  Equation  (2.84)  which  gives  the 

scattered  electric  Field  E  in  the  Fraunhofer  zone  of  the  wall. 

s 

However,  for  purposes  of  comparison,  I|  and  I2  will  be  left  in 
integral  form. 


A- 7 


It  is  obvious  that  the  "Fraunhofer"  integrals  I.J  and  are 
obtained  from  the  corresponding  "Fresnel"  integrals  1^  and  I2  by 

suppressing  the  complex  exponentials  in  the  integrals  of  I.  and  I0 

2  2  L  l 

which  depend  upon  2  and  n  ,  respectively.  The  Fraunhofer  approxi¬ 
mation  is  thus  seen  to  be  based  upon  the  assumption  that  the 
following  inequalities  hold  true  for  all  receiver  points: 


A  =  kh 


(A. 30) 


_  kL^  sin2  (8-h|>)  +  sin2  (y-6) 


(A. 31) 


The  validity  of  the  Fraunhofer  approximation  as  summarized  in  the 
inequalities  (A. 30)  and  (A. 31)  will  now  be  examined. 

In  practice,  the  inequalitily  (A  31)  is  generally  satisfied, 
since  our  computer  programs  automatically  segment  walls  into  sub¬ 
sections  no  more  than  50  feet  in  length.  The  scattered  field  at 
the  receiver  produced  by  each  subsection  is  calculated  using 
either  Equation  (2.84)  or  (A. 24).  These  partial  fields  are  then 
summed  to  give  the  total  field  at  the  receiver  due  to  the  entire 
wall.  This  segmenting  procedure  is  necessary  to  ensure  the 
validity  of  our  assumption  that  f  ($),  the  horizontal  antenna 
pattern  of  the  localizer,  does  not  vary  significantly  as  n  ranges 
from  -L/2  <  n  <  +L/2  and  that  f($)  can  be  approximated  very 
accurately  by  f(ifi)  for  all  points  on  the  wall.  To  get  some  idea 
of  the  magnitude  of  the  parameter  B  appearing  in  (A. 31),  consider 
the  following  set  of  typical  parameters: 

D  .  =  4000  ft 
Pi 

D  -  =  4000  ft 
P2 

L  =  50  ft 

k  =  2r./\  =  .706  (ft)-1 

•The  value  of  k  listed  above  corresponds  to  a  localizer  wavelength 
of  8.3  ft.  Setting  (0+<K'  and  (y-6)  equal  to  90°  (worst  case),  we 
find  that  B  =  .11.  Consequently,  the  approximation  B  <<  1  inherent 
in  the  Fraunhofer  approximation  is  reasonably  good,  and,  conse- 


quently,  differences  between  the  integrals  I2  (A. 26)  and  (A. 29) 
will  generally  be  negligible  in  most  cases  of  interest. 

On  the  other  hand,  the  inequality  (A. 30)  involving  the  para¬ 
meter  A  is  often  violated  for  tall  walls.  Consider,  for  example, 
the  case  of  a  wall  ICO  feet  high  {h=100  ft) .  Using  the  same 
values  for  O  Dd2'  an^  ^  that  were  used  in  the  previous  simple 
calculation  of  the  parameter  B,  we  find  that  A  =  1.77.  Conse¬ 
quently,  since  the  Fraunhofer  approximation  is  only  valid  if  A<<1, 
we  can  expect  to  see  significant  differences  between  the  values  of 
1^  and  and,  hence,  between  the  values  of  the  scattered  fields 
predicted  by  the  Fraunhofer  and  Fresnel  models  when  very  tall 
scatterers  are  involved. 


Specifically,  we  can  certainly  expect  to  see  phase  differences 
between  the  scattered  fields  predicted  by  the  two  models,  since  I, 

A 

is  complex  and  is  real.  Furthermore,  the  intensities  of  the 
scattered  fields  predicted  by  the  two  models  will  differ,  since 
the  magnitude  of  T^(|l^|)  will  generally  be  less  than  the  magnitude 
of  I*(|ljJ).  To  demonstrate  this  fact,  we  first  define  the  fol¬ 
lowing  two  functions: 


g(z) 

In  terms  of  f  and  g,  I 


kz 

2 


(A. 32) 


(A. 33) 


and  can  be  written  as  follows: 


Z1  " 


n 

7 


f  (z )  elgU)  dz 


(A. 34) 


f(z)  dz 


Clearly,  |l^J  must  satisfy  the  following  inequality: 


|f  (z) |  dz 


(A. 35) 


(A. 36) 


A-S 


For  the  small  angles  of  elevation  which  we  have  assumed 
(H/Dd1  <<1,  z2/°p2  <<:  the  function  f(z)  will  generally  be 

non-negative  for  all  z  in  the  range  o  <  z  <  h  even  for  relatively 
large  values  of  h.  That  is  to  say.  in  most  cases  of  interest, 
|f{z)j  =  f(z)  (o<z<h) .  Under  this  assumption,  the  inequality 
(A -36)  becomes 

h 

1 12  i  1  f  f.z)dz  =  =  il{!  {A. 37) 

As  a  consequence  of  the  inequality  (A. 37)  we  can  conclude  that,  in 
most  cases  of  interest,  the  intensity  of  the  scattered  field  pre¬ 
dicted  by  the  Fresnel  model  is  generally  less  than  the  intensity 
predicted  by  the  Fraunhofer  model . 

Figures  A.l,  A. 2,  and  A. 3  show  comparative  plots  of  DDM  in 
microamps,  as  predicted  by  the  Fraunhofer  and  Fresnel  scattering 
models.  The  scatterars  involved  are  flat,  vertical  walls  SO  feet 
wide  and  25,  50,  and  100  feet  in  height  respectively.  The  walls 
are  oriented  parallel  to  the  centerline  of  the  runway.  The  air¬ 
craft  is  assumed  to  be  flying  at  a  constant  altitude  of  50  feet 
down  the  centerline  of  the  runway,  so  that  in  the  absence  cf 
scatterers  the  DDM  would  be  identically  0.  The  abscissa  in  each 
figure  represents  distance  from  the  localizer  in  feet.  In  each 
graph,  the  dashed  line  represents  the  predictions  of  the  Fresnel 
model,  while  the  solid  line  represents  the  predictions  of  the 
Fraunhofer  model. 

Referring  to  Fig.  A.l,  we  note  that,  for  a  wall  only  25  feet 
high,  there  is  very  little  difference  between  ine  predictions  of 
the  two  models.  This  was  to  be  expected,  since  for  small  heights, 
the  parameter  A  defined  in  Equaticr.  t A . 3 C i  is  lih-sly  ::o  be  small 
compared  with  unity,  and,  consequently,  the  predictions  of  the 
Fresnel  and  Fraunhofer  models  should  be  nearly  identical. 

Referring  to  Figure  A. 2,  we  see  that  there  are  slight  differ¬ 
ences  in  both  phase  and  amplitude  between  the  predictions  of  the 
two  models  for  a  wall  50  feet  in  height. 


Finally,  in  Figure  A. 3,  v?e  note  major  differences  in  both 
amplitude  and  phase  between  the  two  models  for  a  wall  100  feet  in 
height.  The  Fraunhofer  model  predicts  DDM  magnitudes  which  are 
nearly  twice  the  correspond  ding  values  predicted  by  the  Fresnel 
model  over  portions  of  the  flight  path.  Differences  in  the  posi¬ 
tions  of  DDM  maxima  of  up  to  two  hundred  feet  can  also  be  observed. 


FRESNEL 


rosnoJ  and  Fraunhofer  Approximations  Compared,  25-Foot  High  Wall 


Figure  A. 2.  FroaneX  and  Fraunhofer  Approximations  Compared#  50-Foot  High  Wall 


Figure  A. 3.  Fresnel  and  Fraunhofer  Approximations  Compared*  100-Foot  High  Wall 


APPENDIX  B 

SCATTERING  FROM  A  VERTICAL  TRIALGLE,  HEM  FORMULATION 


Vfe  use  the  current  distribution  net bod  to  calculate  the 
scattered  field  for  an  infinitely  conductive  triangular  surface . 
The  scattered  electric  far  field  is  (with  ground  reflections  in¬ 
cluded)  : 


iou  e1**2  - 

v  =  £ - a2  x  s2  x  1 

^  «;  *2 


CS-1) 


■where 


I  = 


|c\  1/2  iicRl 

=  8|y|  yfine^sialaOsinCaQ 


.3-2) 


where  A  =  Jc{siny-sinS) 

a  =  kz2/*2 

a  *  ^A^pl 

H,  is  the  antenna  height.  For  a  right  triangle  with  the 
oas^  on  the  ground  and  the  vertex  to  the  left  the  integration, 
dC  is  frees  0  to  h“/3  +  a/2,  and  dn  is  free  -S/2  to  -*-3/2.  The 
following  result  is  obtained  {Figure  3.1} z* 

,ik{P.i+H2)  _  . 

2s 


h  -  ^BEb£l,)S  o 


1  -J— [e1 

|  (a-n)  [ 


2*1 
(m-; a)h/2 


-fR2xk)ccsl 


-  e 


-i(s-n)h/2, 


tincfl 


1  1 

sine  j{A3+0{n-n)h)j 


CA3-C  Ca-n)! 


^r[-= 


>.(»«>  h/Z.Mi 


+  e 


Ca+n) 

-i  (r+n)h/2_ -  I  ± 


sinc^~{A3+o  Cs+n)  h 


sine  I  r{A3-C  (u+n) 


-  “ 
h) 


13.3) 


*Note  that  the  angles  y  and  8  defined  in  Figure  3-1  and  used  bare 
for  convenience  correspond  to  the  angles  ~/2- (~+v)  and  ~/2- (y— 8 ) 
respectively,  used  in  Section  2.5  and  2.6  of  this  report. 
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(3.5) 


sinl3t-r.)h^  sin  (n-rn)h- 
2(3B-al  2te*s) 


1A3/2  iASC 


Khere 


*1  =  A  *  c(=-n)| 
T2  =  A  -  c(a-o)| 
T3  =  A  »  =(n-n  )| 
*  A  OfjfrSSjg- 

C  *  b!/b  -  1/2 


(3.5) 


-nr 


rtex  to  the  left  c=I. 


?cr  the  ease  of  the  vertex  to  the  right  the  integration  &Z 
Is  frcBE  h_  to  -hrt/B  +  fa/2  ass  dn  is  frees  -B/2  to  EC-h^/h  +  1/2)  - 
is  ci^2  by  IB. 4}  ,  •its 
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1  e- 


fct-rJh/2 
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[«b>  v  ^ 

sio(a-o)hL  +  sintatajhjl  T  -iABZ  .-iAS/2  1 
2fe*s)  2(rrr.)  JiA^"  ~  J 


(3.7) 


where  I.  »  T 
So1 
of 


,,  and  Tj  are  gives  by  E.yoat: 
Sgatics  (3.7)  (with  s=-l|  Is  just 
a  (3.5)  *-Itb  —I- 


ifi  G*—ln 

the  ssgls:  conjugate 


These  are  the  new  closed -form  solutions  for  scattering  frcs 
triangular  shapes  that  are  to  be  incorporated  into  the  computer 
program. 


APPENDIX  C 

MULTIPLE  SCATTERING  FROM  VERTICAL  RECTANGULAR  WALLS,  NEW  FORMULATION 

In  this  appendix  we  treat  the  multiple  scattering  of  electro¬ 
magnetic  waves  from  a  set  of  two  vertical  rectangular  walls  where 

infinite  conductivity  is  assumed.  In  the  double  reflection  study 
4 

given  by  I.B.M.  only  the  reflected  field  was  used  in  their  Equa¬ 
tion  (2.3).  For  reflection  from  a  reflecting  rectangular  wall, 
the  current  distribution  method  yields  the  same  result  as  given  by 
I.B.M.  if  the  total  field  is  used  instead  of  the  reflected  field 
and  if  the  second  line  integral  in  the  right  side  of  their  Equation 
(2.3)  is  included.  In  the  following  calculation  the  current  dis¬ 
tribution  method  is  used  to  obtain  the  new  equations  for  double 
reflection. 

The  electric  field  at  a  point  on  the  second  wall  is  given  by 
(Fig.  C.l)* 


En  —  i 


kLlhiEpf (») 
"R1R2 


exp|ik(Ri+Dp20+n2sinY2)j 


(C.l) 


where  sinc(x)  =  sinx/x  and  where  k  =  2tt/A. 

At  the  second  wall  the  incident  magnetic  field  is  given  by 

Hp  5  /iTu  Ep(-Jc)  (C.?.) 

and  the  surface  current  on  the  wall  is  given  oy 


♦Note  that  the  angles  yi  and  £5.  defined  in  Figure  C.l  and  used  here 
for  convenience  correspond  tothe  angles  tt/2  -  (0+^)  and 
tt/2  -  (y~9)»  respectively,  used  in  Sections  2.:  and  2.6. 


C-l 


m . . 


where  (with  ground  reflection  included) 


The  scattered  electric  field  at  the  receiver  from  the  second  wall 
is  then  given  by 


-*•  *  a  k2  ,  exp  ( ik  [Rl+Dp2a+ft3]  ) 

E„  «  -  iR-sfik  -rf/fi^hjEof  (** - - - ““ - *-cos&2 

s  3  A  z  x  °  RlI>o20R3 


where 


fL2  1  ftl  1 

cosBiosinc^J^^ts*  :Y2~sin82)J  sinc^— k(sinYX-sinGXQ)j 


sin (AC) 
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C-h. 
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2  cos (A+B)x  C+h2  cos(A-B)x 
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dx-  f 
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dx 


A  =  kx^/R^  * 


_  khx  =  khx 

Dp2  Dp2o' 


=  DP20H 
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and  where 


/*sin  ax 
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(ax)- 


(C. 


(C.l 


(C .12] 


/cos  ax  (ax)z  (ax)s 

__d*  *  log  (ax)  -  ^r-  ♦ 

Combinations  of  reflections  from  walls  with  the  lower  edge 
flush  against  the  ground  plane  can  be  used  to  compute  the  reflec¬ 
tion  from  walls  whose  lower  edges  are  at  a  height  above  the  ground 
plane.  The  following  figures  show  schematically  such  a  procedure. 
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Figure  C.3.  Elevated  Structures 

SYMBOLS 

from  antenna  to  reference  point  on  ground  plane  at  1st  wall, 
from  foot  of  antenna  to  reference  point  at  1st  wall, 
f  (ifi )  antenna  gain  factor, 
h^  height  of  1st  wall, 
width  of  1st  wall, 
angle  of  incidence  at  1st  wall. 

look  angle  from  reference  point  on  1st  wall  to  element  of 
integra.  on  on  2nd  wall. 


6lo 

A 

*2 


Dp20 

D 

P2 


n 


■3 


look  angle  from  reference  point  on  1st  wall  to  reference 
point  on  2nd  wall. 

normal  vector  at  1st  wall. 

vector  from  reference  point  on  1st  wall  to  element  of  inte¬ 
gration  on  2nd  wall. 

reference  point  on  1st  wall  to  reference  point  on  2nd  wall. 

reference  point  on  1st  wall  to  foot  of  element  of  integration 
on  2nd  wall. 

width  of  2nd  wall. 

height  of  2nd  wall. 

normal  vector  at  2nd  wall. 

unit  vector  normal  to  ground  plane. 

2u  divided  by  wavelength. 

unit  vector  along  horizontal  direction  tangent  to  wall  surface 
at  2nd  wall. 

vector  from  reference  point  on  2nd  wall  to  the  element  of 
integration 

angle  of  incidence  at  2nd  wall, 
look  angle  from  2nd  wall  to  receiver. 

vector  from  reference  point  on  2nd  wall  to  aircraft  receiver. 


C-6 


APPENDIX  D 

DERIVATION  OF  GAIN  VECTOR  FOR  SMALL  CIRCULAR  LODP  RECEIVING  ANTENNA 

We  consider  a  linearly  polarized  plane  wave  given  by 
E  =  EQe  i(k*r-wt)  incident  on  a  perfectly  conducting  circular  loop 
of  radius  a  which  is  connected  to  a  receiver  input  impedance  Z 

L 

(Fig.  D.l).  The  curren.t  generated  in  the  load  impedance  2  can  be 

J  JLi 

written  as 

I  =  V/(Zl  +  Za) »  (D.l) 

where  Z  is  the  antenna  radiation  impedance  and  V  is  the  net 

cl 

electromotive  force  (E.M.F.  developed  in  the  loop  by  the  incident 
field).  The  E.M.F.  is  given  by  the  following  development  of  the 
basic  line  integral  of  E  around  the  loop: 

V  =^*E  •  dl 

■  ff  dA  n  •  ^xE 

.  -*•  -*■ 

=  iEQ  •  (nxk)  ffz  dr  do  elic  r  (D.2) 

Here  the  vector  r  lies  in  the  plane  of  the  loop  and  is  perpendicu¬ 
lar  to  the  unit  vector  n.  We  may  replace  the  scalar  product  in 
the  exponential  by 

k  •  r  =  k  r  sinS  cos* 

Q 

with  the  result  that  the  integral  reduces  to 

r3  j  ja  _2kr  sinS  cosO 

f  r  dr  J  dO  e 

=2“  f r  dr  J (krsinS) 

•f  o 

2 

=  2-a  (ka  sinS)/(ka  sin5) 

=  »a2  (D.3) 

The  latter  approximation  is  valid  whenever  a<<X/2-.  Using  this 
re;. alt  we  find  that  the  current  delivered  to  the  receiver  can  be 
expressed  as 


D-l 


I  = 


ita* (nxk) *E0 


ZL  +  za 


(D.4 


Thus  the  gain-polarization  vector  of  Equation  (D.l)  takes  the  form 


g(k)  =  can  x  k 


Figure  D.l. 


Illustration  of  Circular  Loop 
Receiving  Antenna  Problem 


(D.5 


This  vector  selects  the  component  of  E  normal  to  the  vertical 
olane  of  incidence  and  aoolies  a  factor  of  sinS. 


APPENDIX  £ 

PROGRAMMING  ASPECTS  OF  THE  ILS  MODEL 


I.  INTRODUCTION 

The  complete  simulation  of  apparent  CDI  for  a  large  group  of 
ground  structures  is  accomplished  by  using  a  sequence  of  programs. 
First,  the  DATGEN  set  of  routines,  written  for  a  time-snaring 
PDP-10  computer,  are  employed  to  set  up  descriptions  of  antenna 
patterns,  receiver  trajectories,  scatterer  geometry,  etc.  in  the 
correct  formats  for  input  into  the  ILS-O-E  program. 

The  ILS-O-E  program  is  the  basic  simulation  program,  which 
can  be  used  to  calculate  static  or  dynamic  CDI  along  specified 
receiver  trajectories  for  a  variety  of  antenna  and  scatterer  con¬ 
ditions.  The  cutout  of  this  program  is  usually  a  formatted  listing 
on  tape  of  CDI  versus  distance. 

If  the  static  CDI  has  been  generated,  the  dynamic  CDI  may  be 
obtained  with  program  ILSDYN  without  further  recourse  to  ILS-O-E. 
The  CDI  data  cure  then  usually  input  to  the  ILSPLOT  program  for 
final  presentation  in  graphical  fora.  The  ILS-O-E  program  is  a 
complete  operational  unit  in  itself,  and  its  use  is  detailed  below 
in  terms  of  the  input  data  specifications.  However,  we  do  note 
that  since  the  ILS  analysis  program  and  the  ILS  computer  program 
are  presently  under  development,  the  ILS-O-E  computer  program  is 
considered  an  interim  and  preliminary,  though  fullv-operational 
model . 


II.  DESCRIPTION  OF  DATA  GENERATION  PROGRAMS  FOR  ILS-O-E 

A  series  of  programs  have  been  written  on  the  PDP-10  to 
generate  input  data  for  the  ILS  program.  The  first  of  these  pro¬ 
grams  is  STRUGN.  This  is  a  Structure  Generator.  It  takes  section 
descriptions  in  free  field  format  and  produces  a  list  in  the  proper 
format  for  input  to  ILS.  The  sections  may  be  described  singly  or 
as  polygons  or  as  sections  of  circles.  In  the  latter  case,  the 
arcs  will  be  approximated  by  a  specified  number  of  chords.  Any 
number  of  sets  of  sections  may  be  given.  Each  set  is  considered 
a  structure  and  is  delimited  and  numbered  in  the  output  list. 


The  output  of  STRUGN  {or  a  combination  of  the  output  of 
several  runs)  is  used  as  the  input  to  FLDGEH.  This  is  the  Field 
Generation  Program.  FLDGEH  reads  in  each  structure  in  an  input 
list.  It  then  types  out  the  number  of  the  structure  and  accepts 
from  the  teletype,  input,  specifying  the  locations  and  orientations 
of  the  structure.  Each  structure  nay  appear  any  number  of  tines 
on  the  field  or  nay  be  emitted  completely . 

This  field  description  is  used  as  the  input  to  LOCAT .  LOCAT 
takes  the  field  and  arbitrary  antenna  location  and  orientation  and 
determines  where  the  sections  are  from  the  antenna.  For  each  sec¬ 
tion,  LOCAT  examines  the  relative  angle  of  the  line  of  sight  and 
the  normal  to  the  section  to  decide  if  the  section  face  will  receive 
any  incident  field  from  the  antenna.  It  then  outputs  a  list  of  the 
sections  making  up  the  field.  If  a  section  is  not  illuminated,  a 
flag  is  set  in  the  output  to  indicate  this.  As  LOCAT  can  be  run 
with  the  same  field  data  and  different  antenna  locations,  prepar¬ 
ing  data  for  an  airport  with  multiple  antenna  locations  is  simpli¬ 
fied. 

For  simple  cases,  the  output  of  LOCAT  can  be  used  as  is,  but 
for  involved  airports,  there  is  a  problem  of  shadowing.  A  structure 
section  may  be  facing  the  antenna  location  but  not  have  an  incident 
field,  because  another  structure  blocks  it.  To  determine  if  this 
is  the  case,  a  series  of  programs  (SORT  1,  SORT  2,  SORT  3,  PRUN) , 
hav.e  been  written.  These  programs  take  the  output  from  LOCAT  and 
find  out  which,  if  any,  sections  are  shadowed  and  create  a  new  list 
with  flags  to  indicate  the  hidden  sections.  The  final  program  in 
the  series  eliminates  those  that  are  flagged  as  shadowed,  as  well 
as  those  flagged  by  LOCAT.  This  results  not  only  in  an  increase 
in  accuracy  of  the  ILS,  but,  as  the  hidden  sections  do  not  need  to 
have  their  scattered  fields  computed ,  also  in  a  substantial  savings 
in  time. 

To  simplify  verification  of  the  input  data  used  in  the  pro¬ 
grams,  a  plotting  program  was  written.  This  takes  the  output  from 
FLDGEH,  STRUGK,  LOCAT,  PRUN  or  input  cards  for  the  7094  and  draws 
the  sections  in  their  proper  locations  and  orientations.  The  scale 


E-2 


is  variable  for  convenience  in  comparing  to  blueprints.,  f iocrplans 
etc.,  and  hidden  sections  cay  be  indicated  by  dotted  lines,  if 
needed . 


III.  IKPJJT  DATA  FOR  ILS-C-E 


The  first  input  card  is  the  carp -station  cede  card.  It  has 
the  following  f oread: 

Col.  Symbol  Contents 


1-2  Mode 


3-4  2DFC 


11-20  FRQ 

21-30  TEETGS 

31-40  XTE 

41-50  Slope 

51-60  ZA  (1)^ 

61-70  ZA  (2)  l 

71-80  ZA  (3) J 


1  (Y-  Ring) 

2  {8  Loop) 

3  (Waveguide  localizer) 

4  (i*ot  used  at  present) 

5  (Single  censured  localizer  pattern) 

6  (Measured  capture  effect  localizer 
pattern) 

C  ( Fraunhofer  approximation} 

1  (Internal  switch  for  Fraunhofer/ 
Fresnel) 

Frequency,  MH2 

Glide  slope  angle.  Beg. 

Distance  free  ice.  antenna  to  threshold 
Inclination  of  runway  at  threshold,  deg. 

Antenna  decent  heights,  ft. 


The  code  selects  the  type  of  antenna  pattern  used  for  the 
simulation .  For  a  single  frequency  system,  a  theorecticai  pattern 
(V-ning ,  8-Lcop,  Waveguide)  or  a  measured  pattern  nay  be  used. 

Two  patterns  are  used  for  a  dual  frequency  sy, cem,  such  as  the 
Alford  capture  effect  system.  The  clearance  pattern  is  chosen 
first.  This  is  done  by  using  a  negative  value  for  the  node  number 
If  a  measured  pattern  is  desired  for  the  clearance,  a  node  of  6  is 
used.  After  the  clearance  pattern  is  set,  the  course  pattern  is 
chosen.  The  IDFC  switch  is  used  to  force  the  simulation  to  use 
the  Fraunhofer  approximation  if  desired;  otherwise,  the  Fresnel 


■ertical 


Los  will  be  used  for  those  surfaces  vita  large  enough 
to  it- 


If  an  aatenna  catte: 
Is  used: 

Col. 

1-10 

11-20 

21-30 


as  to  re  reao 


f: 


G  rcrsaf 


Ar.de  of  the  reasurenent  ia  degrees 


as  sea saras 


Csaasite  sac 


usaosite  earner  patters  as  seasare 


ihe  assies  east  be  between  plus  ass  sinus  ISO  degrees  and  in 
ascend  ins  order.  A  mximr  of  fifty  aeasaresaents  is  allowed.  If 
less  than  fifty  are  required,  the  last  seasaressat  should  be 
followed  by  an  angle  is  excess  of  4GS  degrees. 


The  next  card  is  the 


Col. 

1-10 

11-20 

21-30 

31-40 

41-50 

51-60 


Svnhol 


surse  -.lath  ears  in  tse 
Contents 


Llowms  rc 


aSA  (1)  Distance  along  runway  centerline  to 
course  array 

m  C2)  Distance  along  runway  center 1 ire  to 

clearance  array 

•Sot  Uses! 

CW  Course  width 

fN'ot  Osedl 

CDS  Relative  clearance  signal  strencth 


If  a  course  width  of  greater  than  three  degrees  is  read  in, 
this  width  is  us*-d-  If  a  width  of  less  than  three  degrees  is  in¬ 
put,  the  width  wxli  be  adjusted  to  FAA  specifications  according  tc 
the  distance  to  the  threshold,  as  giver,  on  the  node  card. 


The  recei 

pdth 

card  follows  in  the  : 

Col. 

Symbol 

Contents 

1-10 

XMBi 

Starting  distar.o 

11-20 

XXAX 

Ending  distance  : 

21-30 

DXR 

Soacing  between  : 

31-40 

PEIR 

Angle  of  apprcaci 

41-50 

?SIR 

Glide  path  ancle 

51-60  (Mast  be  left  blank) 

Si -70  ZuP  Height  above  glide  path 

71-30  (Must  be  left  blank) 

If  sore  --has  50G  receiver  points  are  requested ,  the  program 
vill  increase  DIB  until  rot  sore  than  500  will  be  needed.  The 
nusber  of  points  is  equal  to  XMAX  nines  IMIS  divided  by  DXA.  If 
a  clearance  ran  is  desired,  the  format  of  the  card  is  as  follows: 


Col . 

Symbol 

Contents 

1-10 

xmo 

Initial  angle 

11-20 

XAXD 

Final  angle 

21-30 

dxe 

Angular  increment 

31-50 

(Must  be 

left  blank) 

51-60 

BO 

Radius  of  orbit  in  feet 

61-70 

£IIF 

Height  of  orbit 

74 

XCF 

1  (Switch  to  indicate  orbit  run) 

Tbe  next 

card  gives 

the  velocity  of  the  aircraft  in  feet  per 

second  (colmms  11-21) .  As  a  zero  velocity  will  cause  erroneous 
results,  the  program  sets  a  sinistra  velocity  of  one  foot  per  second. 

The  next  series  of  cards  descr*be  tbe  scattering  surfaces. 


Each  surface 

is  a  single 

vertical  or  tilted  flat  wall  of  infinite 

conduct  ivity . 

Col. 

Symbol 

Contents 

1-2 

ID 

1 

3-8 

XW  (1) 

X-Coordinate  of  center  of  base 

9-14 

X w  (2) 

Y -Coord inate  of  center  base 

15-20 

XW  (3) 

Zero 

21-24 

Rone 

(Blank) 

25 

C 

Zero 

26-30 

ALPHA 

Angle  between  base  line  and  X-Axis 

31-35 

DELTA 

Zero 

36-50 

(Must  be  3 

.eft  blank) 

51-60 

MW 

Width  of  wall 

61-70 

EW 

Height  of  wall  free,  top  to  bottom 

71-30 

EL 

Height  of  base  above  ground 

Any  number  of  surfaces  say  be  used,  limited  only  by  available 
computer  time-  Complex  surfaces  nay  be  approximated  by  using  sec¬ 
tions  as  chords  of  arcs,  etc.  If  sections  exceeding  50  feet  in 
width  are  used,  they  will  be  broken  up  into  sections  of  less  than 
50-foot  width  by  the  program.  Thus,  no  tine  will  be  saved  by  using 
approximating  pieces  wider  than  50  feet,  and  accuracy  will  suffer. 

If  an  ID  of  -1  is  used,  the  scattered  field  is  subtracted  from  the 
total  field.  Thus,  if  a  wall  is  represented  by  one  or  sore  s c bit¬ 
terer  s,  a  hole  in  that  surface  can  be  represented  by  a  series  of 
scatterers  with  an  ID  of  -1  - 

After  a  series  of  scatterers  is  finished,  a  final  ID  is  in¬ 
serted.  An  3D  of  10  will  cause  the  direct  ray  to  be  added  to  the 
field  summation  and  a  CDI  tape  tc  be  generated  for  the  plot  pro¬ 
gram.  Then  a  new  series  of  scatterers  say  be  input.  An  ID  of 
11-15  will  cause  a  new  receiver  path  to  be  read  in,  allowing  multi¬ 
ple  runs  in  one  batch.  If  the  ID  is  less  than  14,  the  plot  program 
will  combine  the  runs  into  one  graph.  This  allows  flight  paths 
containing  core  than  500  points  to  be-  sade.  However,  this  will 
require  th^t  all  the  scatterers  be  resubmitted  for  each  portion  of 
the  flight  path.  On  any  one  section  of  the  flight,  the  receiver 
point  spacing  must  be  constant,  but  between  parts  of  the  flight 
different  spacings  say  be  used.  This  will  allow  substantial 
savings  in  computer  time,  as  those  areas  with  slow  changes  in  CDI 
(such  as  near  the  outer  marker  and  beyond)  will  not  require  the 
density  of  points  that  areas  with  rapid  changes  in  CDI  (such  as 
near  the  threshold)  will. 

At 'the  end  of  a  series  of  runs  on  a  flight  path,  an  ID  of  20 
say  be  used.  This  will  cause  the  program  to  read  in  a  new  mode 
card  allowing  the  program  to  start  over  with  a  new  antenna,  flight 
path,  and  scattorer  set. 

An  ID  of  0  will  cause  the  direct  ray  to  be  calculated,  thus 
permitting  CDI  measurements  of  the  undisturbed  field.  Khen  this 
ID  is  used,  not  only  the  CDI,  but  the  magnitudes  of  the  carrier 
and  sideband  fields  are  output.  This  allows  verification  of 
antenna  patterns  and  comparison  of  thecrectical  and  measured  results. 


E-6 


To  terminate  the  program,  an  end-of-file  card  is  inserted. 
This  will  cause  the  program  to  output  the  CDI  if  an  unfinished  run 
is  in  progress  and  then  stop. 
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